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Abstract 

Sludyiufi; (iu* dynamic response of a car suspension system and to find the ojitimal response in 
vai itiiis kinds of t'xcilations is one of the major challenges that the automobile iialuslry is facing. ‘l’in' 
designer has to strike a balama* Ix'tween a vast number of obje<'tives. many of thc'in coiiilicling to each 
Miller. subj«‘cl to si'vi'fal constraints and obtain an optimal solution which does not ttverlook any of the 

I il ijeci ives. 

rh(' primtiry function of ti stispension system in a car is to isolate the rotul excitations experienced 
Ity till' wheels from being translern'd to the passtuigers. In the present study we tire considering an 
iiuiependent suspimsion systtun and we are assuming the part of the car abovi' the chtisis as a rigid 
body, 'riie suspension system is modelled as two and threi' dimensional dynamie tnixiels. Tlu' car body 
is iisnally snpporti'd by a snsjtension spring and a damper at eaeh wheel. The tyn* of the wheel can be 
;dso .assumed to have some stiffness in the vertical direction. The dynamic behavionr of a ear snsiten- 
■•'ion system is di'scribed by nonliueai' differential equations of motion found from the frtw^ bexly ditigram 
aitproach. The design vtiriables ar(^ the stiffness and the damping parameters. Other ptirtimetcM's. which 
.di'ecl tlie dynamics of the car motion, like wIkhiI to wheel distance, front and rear unsprung masses etc. 
.lie taken as constants. Different tyi)es of shocks, i.e., road excitations are ajjplied to th<i whet'l such as 

I I lose caiisetl by the wheels striking a bump, whetils striking a series of bumps, and that causeil by the 
wiii'ers falling into a pot-hole et<'. 

A real world problem such as this poses many problems to the d^igner. Mast of the conventional 
ie( lmi(|nes are not suitable for this problem, as apriori knowledge of the car dynamics is recpiirtxl. Also. 
I li(!si! processes are very time consuming and computationally expensive. Besides, there is a high prof)- 
ahility that the results obtained are local in nature. The use of genetic algorithms overcomes the above 
mentioned diffietdties and so the present study entails the use of genetic algorithms. 

¥ 

'lb .sttirt with a single objective such as maximum bouncing, pitching or rolling mnplitude of the 
s|)nmg mass is minimised subject to satisfying a number of constraints. The objective is expressed in 
t erms of ’transmissibility’ which is the ratio of maximum dynamic amplitude to the maximum excitation 
amplitude. Also, the objective function is considered in two different ways. Firstly, for optimization 
in-ocess we consider the time interval after the car has cleared the bump and in the second case we 
consider the whole period of travel including that on the bump> The constraints mainly arise from 
(■(imibrt considerations such as limiting the maximum vertical jerk experienced by the sprung mass, 
frotil. and rear natural frequencies, and others. Even for a single constraint, restricting the lufvximiim 



\ (n i.i(vil jcirk. the feasible search space is found to be complex and multi-modal. 


A real-coded genetic algorithm can successfully find near optimal solutions for different objectives. 
Since? tli(! sf'arch space is real the precision is more as compared to a binary coded genetic algorithm.The 
snlni ioiis obtained are better than a design that was being used in an automobile industry. 

How(!ver. the solutions obtained by single objective optimization are not fully satisfying for practical 
ap])lications. Dynamic system design can be formulated as a multi-objective problem. In the case of a 
ca r susi)eusiou design, the designer is faced with conflicting objectives such as minimization of maximum 
vcrlical displacement and minimization of maximum vertical acceleration. It has been observed that 
1 lie smaller maximum vertical displacement leads to larger value of maximum vertical acceleration. So 
it is desirable to obtain a set of optimal solutions all equally good, firom which the designer can choose 
a particular solution according to his requirement. 

Non dominated sorting genetic algorithm-II (NSGA-II) can efficiently tackle the above multi-objective 
|)rul)l<-m. This algorithm gives us a set of optimal solutions which takes care of both the conflicting 
nl)j('ctives in a .single run. 

Under steady state excitations, the damper can be modelled as a linear system. However as the 
si lock loading comes into picture the nonlinear property of the dampers come into play. So a non-linear 
i iihic damping term was considered over and above the usual viscous damping term when the car was 
Milijected to a rounded step displacement function. 
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Chapter 1 


Introduction 

1.1 Introduction 

I Jcsif^n (if (lyiijunic sy.sl.<'ins is usually Inised on intuition and experience. Although industrial couipaiiics 
lia.v(' slartcKl to switch from experimental studies to computational methods l)a.scd on the multi-body 
sysicm approach, (h^slgn of mechanical sy.stems is still largely carried otit by intuitive* chaug<'S of the 
di'sigii vfiriabh's. This is due to the fact that the multi-body .system approacth is well develoiu'd I'or 
iiialyziiig the dynamic Ixdiaviour of mechanical .systems iii a variety of disciplines, btit there is a lack of 

I I Mils for a systi'iuatic design via optimisation. The optimisation methods an; not so popiihir with the 
ilesign ('ugiiKstrs in tlu' industry mainly because most of them ai’e not having wide range of us<*. Th<;y 
iisiitilly work v(ny well for a certain kind of problems and fail miserably in tint other ciuscts. Sr>cou<lly. 
most of them nsjuire a certain amonut of apriori kiiowhtdge abotit the problem whit^h is not always 
• ivailalile to lh<( designers. Thirdly, these methods requin; mathematical modelling of the nsil wot Id 
jii iibh'nis whi<;h is not always <;fisy. Also, these methods are very sensitive to th<* initial .solution chosen 
.iiid h(uic<f retiuint a number of runs to obtain a global optimal solution. This mak(!S the.so methods 
very time constiming and computationally expensive. However, genetic algorithms, which work on the 
iiMt,ur(('s principle of the 'nunnval of the fittest', can overcome all the above shortcomings of (tonventiomil 
I i|)t imizatiou im?thods. In the present study the dynamics of a car suspension have been o[)timi5!ed using 
genetic algoi'ithms. 

Modelling and simulation in the field of vehicle dynautnics is a complex topic. For the optimization 
of the vehicle’s dynamic behaviour, not only model accuracy but also an efficient implementation of 
suitable approaches is of great importance. Pra<ctical work on this topic has resulted in integrated 
( ffiatmeut of the mathematical modelling combined with adapted integration schemes. This approach 
also seems to be very well suited for applications in optimization. 

Suspension systems are used to isolate the road excitations firom being tramsmitted directly to the 



p.issi'iijiitTs. S<*«-oi»Hy thoy nmsl do it without impairing tlio staiulity. stporing or general liandliim t 
ii ii's of tlu' voiiiclc!. A modorn vohhde suspension system is n complex oscillalory sysleni with a niunhei 
mI I'lasto-daniping elements, the characteristics of which are in most ca.ses of nou-linear <'haiai-lfr. Dili' 
111 I III' car’s motion on the road and also to the eftects of the power-train and the rot at ing elements, i he 
M'liiclc jierfonns stochastic, secondary, .spatial oscillatory motion that affects tin* stability and coinfoii 
|)aramet<'rs. and may leail to fatigue of vital parts of the car. which in turn may threaten traflic .safety. 
The negative ('tt'erts of the secondary ascillatory motion can be redu<-ed by way of optimization ofehar- 
aci eristics of the ('lasto-damjiing elements of the vehicle. 

.\lt hough ililferent types of suspeihsion systems exist, we consider here the design of an indepeiideiii 
Mispension systi'in [1]. The snspimskm system parameters liki' stifiue.ss of the .springs iind the damping 
coellicients of the dampers are treated as design variables.Tlie optimization problem eonsiders maxi- 
mizing the’ ridi' comfort by minimizing the trausmissibility and acceleration of th(> sjirung mass when 
I he sns]K‘nsion .system is subjected to different kinds of base excitations. 


in 1 lie present study, we ttonsider both the two dimen.sional model and three dimensional model ol'thc 
car. Firstly, we use real-coded geneth; algorithm (RGA) [2], a single objective optimization algorit hm, in 
conjunct ion with mlaptive step-size control Rnnge-Kutta method [3], a numerical inU'gration technique. 
I'lii' liounce motion of the car is considered and the minimization of trHn.smis.sibility is consideri'd to be 
liic ulijcctive. .‘\ (uinstraint limiting the value of nuvximum vertical jerk is u.se<i. The uhtaiueil results 
.tic compared wit h previously (jxisting values of the design variables. 


.Since (rm'-snspension design requires optimizing multiple objectives, many of which are ctaiHicting 
to um? another, it is desirable to have an optimization method which can effectively haudl<> the.se objiH> 

I i\'es sinmIta,neousIy. In the present study, we use the multi-objective genetic algorithm. Non dominatcfl 
sort ing genetic algorithm -II (NSGA-II) [4], for the purpose of optimization. It is u.sed in conjunction 
with a.da]>tive .st<ip size control Runge Kutta method, a numerical integration techni(}ue. 

NSGA-II algorithm considers minimizing the vertical amplitude and minimizing the vertical accel- 
eration as a set of conflicting objectives subjected to vertical jerk and frequency constraints. In case of 
thri'e dimensional model, we also consider other set of objectives like pitching amplitude and pitching 
accideration, rolling amplitude and rolling acceleration. The constraint limiting the maximum allow- 
able vertical jerk experienced by the sprung mass is considered. The front, rear and pitching natural 
frequencies are limited to a lower and an upper bound. Finally in order to dampen the pitching Kiotion 



l.isii r. I Ilf n.Mtmal tVfcjnoiu'.v of thf front susjx'nsion systoni is kopt snmllfr tluui ihaf of ihf icai . Alsu. 
I iif vai iahlf liuuiids am fixed in order to avoid extra soft suspeuskm. Tliis leatls to Ik'IIit liaudlin<; ami 
ruinrol of tli(> ear. Different kinds of road exeitations like a sinusoidal Innnp. a series of Inimps. ami a 
loiimled stej) displai’enient function etc. are considered. To tfike into a<‘comit the iran.sk'iil re-sponse 
I iiaracterislics of a damper, cubic non-linear daminng term is considered in uihiilion to the usual vi>- 
c oiis damping term [()]. The successful application of NSGA-Il suggests that it is an ellicietil imuliod 
lor handling dyuainic .system design i)rol>lems in mechanical euginetiring. The same soliwiire can be 
iiM'd for diffon'iil dt'sign objectives for a (tar stuspension by a small modification in the objectives and 
the fonstraiuls definitions. 


1.2 A brief overview of previous work 


('.X. Spentzas n.scHl Box's method [7] to determine the (dastic and damping characP'ristics of a ve- 
hicle susiiension that ojitimize the ride characteristics of a ctir. A sev(m degrees of frei'dom vehicle 
model was used and springs and dampers with non-linear chai’ticteristics wore considc'red. Objective 
I'nnction considentd vertical and rotational acceterations of the vehicle, permitting the optimization of 
ride eharact.eristics. .-Mso. by adding a term in the objective function related to the dynamic tyre load 
(.trial ion. the handling behaviour of the v(;hicle was simultan(K)usly optimiztMl, Only the rattle space 
constraint was considered. Only one objective wfts considerexl while optimization and other confiictlng 
obiectives wert' ignon'd. The jerk constraints (for comfort ride) and the fretpieucy constraints w(*rt.i not 
considered. Also, the process is computationally expensive and time (xmsumiug. Using Box's method 
lor opt imization of non-linear springs and non-linear dampers, with so complex a model recjuin'd a 
powerful computer and long computation time. Consequently, the use of scwen degnws of frecHloni 
model was not justified from an engineering point of view. So, the use of four degrees of model was sug- 
g((st(!d. while still considering non-linear springs, non-linear dampers and large oscillations of the vehicle. 


Miroslav Domic [8] optimized the characteristka of the elastic-damping elements of a passenger car 
by means of a modified Nelder-Mead method (Rao, 1984). In the work reported in the paper, tlu' author 
modified the method of Nelder-Mead by introducing design-dependent bounds on the oscillatory param- 
et.i'r values. These values were computed on the condition that the object function has minimum, that 
tlu! resonant frequencies of the aforementioned masses are located within the recommended bounds, and 
i.hal, the oscillatory parameters are within the prescribed bounds. The parameters computed provide 
for optimal comfort of the vehicle considered when driving on asphalt roads at high speeds of about 30 
m/s. He applied this method to the plane model of a passenger car in order to minimize simultaneously 
v('i tical oscillations of the front and rear wheels, pltdi, and vertical oscillations of a partially loaded 



• i lii. Ic, Xi) lVi-(n»‘ucy <>i' jerk ctnistraints \v<'ro dioson. FnriluT. .sin<‘(.‘ siiiipk^x scari'li iiictliii(l iisfil. 
I iic uliiaiiu!(l sulntidii wjus ck'jx'iidonl <m tlu' initial solutions. Tlion? is a hijnh proltability llial tlio altin- 
1 it Inn may Intul to ti lortil optimum. 

In anoth(-r papi'r. Domic [7] din'ctcd his imtfmnutc towards a modified vcrsitni of (In' Ihtokc-.lcovi-s 
iiii'i IkiiI. .scokiitu; the optimiztition t»f tlie characteristics of the elasto-danipiu}!; elenn'iits of (■ars from I In* 
l not only of ride properties, hut also of handling. Delimitation of tlie range of the opfimi/ation 
liar.nni'lcrs was achieved by imtorpornting con.straints to the algorithm by means of e.K tenia 1 penalty 
Innctions. Tlu^ vehicU' mothd that was used was obtained by extending the ithissic seven degrees of free- 
dom mod('l (bomute, roll, pitch motion of the sprung uias.s an<l bounce of the four unsprung masses) to a 

0 11 degrei' of freedom one. The three .supplementary degrees of freedom (forwtinl displacemmit . latinal 
lisplacement and yaw angle) permit the handling behaviour to be taken into consideration. But. since 

1 lie matrix of transformation of co-ordinates from the body-fixed co-ordinates system to the groimd- 
lixed co-ordinate system was linearized, only small displacements of the vehicle could be considered ami 
I lie study of extreme* situations either for the ride behaviour or for handling is not possible. 

.Another approach to this prolilem was prtwentwl by Pintado and Benitez [7] who ns<*d the lini'ar 
I li t igramming method. The problem was solvetl in the time domain and also in the frecjtumcy domain, for 
detfiniinistic and random excitations. But their appmach coiwidered only liuoar .springs and dampers 
anil this fact, restricted the practical interest of tlwi method, since mtjst cars now-ti-days n.s(> non-linear 
springs and damper.s. 

Wiminer and Ranth [9] formulated the car suspension .system as a multi-criteria optimization of 
maximizing .safety as well as maximizing the ride comfort. The optimal points were found by scalar- 
izat.iou of the vector optimization problem. Frequently used sccUarization approaches are the weighting 
n il.i'i ia method, the global criterion method, the weighting min-max method, and the hierarchical op- 
limizatioii method, respectively. The remaining single criterion optimization problem was solved by 
standard optimization algorithms. However no aspects related to the jerk or frequency of oscillation of 
I he i)a.sscnger.s were considered. 

Markiue et al. [10] considered a finite element model of a trude and formulated an optimization 
in-oblem with the location of various elements at the suspension system as design variables. The dy- 
namic behaviour of truck was approximated by linear finite element models (both 2-D and 3-D). The 
I I lad surface profile was presented as a random function with known power spectral densit;^. The design 



aii:ilil(‘s (•(»uiprisi> gtHdiH^try as wall as spring and damper j)roperti(\s. LiniUaliuiis were inipuM'd un 
iiiaxiiiniin values of t he relativ<‘ displaeonients of suspensions, dynamic wIkh'I loa«l fti a.xles. aeeeleraiiiui 
Ilf the car. A inulti-iH)int approximation methcxl was used. 

.Also. Dieter Bestlci and P(^ter Eheluird in their paper [1 1] c.onsidered optimizat ion t)f a v('hicl<- siispeii- 
>iiiii system, 'riiey eousidiTod two conflicting objectives, pitching ac<Htlerat km an<l vtu'tical acc(‘ler:itioii 
(11(1 used non linear programming along with a numerical integration teclmiqm* to obtain optimal so- 
liiiiiius. However, they u.s('d the weighted objectives methotl as a ntpresentative of the scJilarizal ion 
jii iiiciple. Although here the hierarchical .sti’ategj' converged to optimal points, the problem of local 
iiiiiiiiua existed. 

.None of tilt' above studies discus.sed the complexity of the feasible search spact? in a car suspension 
design problem. Due to nonlinear interactions of design variables and the retiuiremeut of .solving a std of 
iioiiliiiear ililferent.ial equations governing the car, the search spaxie in the.se problems is usually complex 
.iiid multi-modal. 

Dt'b and saxena ustxl a simple binary codtxl genetic algorithm [1] }us an optimization tool to find 

I he optimal susiiension .system for a number of objectives stich as minimizing bouncing, ijittthing find 

I I illiiig funiilitiuies. Frequency fuul jerk constnunts wore also considered. The mlvantage of using geimtic 
algoi'ithms wfus manifohi. It cfui lie cfusily compounded with a numerical integration technkjut' to obtfiiu 
the optinnil solutions. Secondly, unlike other conventional techniques this algorithm starts with fi set 
of solut ions and does not mislead us to a local optima. The solution obtained is fi global solution and 
i.-^ iudepeiideut of the initial .set of solutions. Thirdly, no derivative information is required to rea(;h t(j 
I he opt ima. Also, it is computationally lesser expensive and lesser time con.suming. 

However solutions obtained with such an approach are not fully satisfying for practical applications, 
since they take into consideration only a single aspect out of a couple of conflicting system requirements 
oi- objectives. The present work deals with the car suspension optimization problem in a multi-objective 
(•(intext. The multi-objective approach seems to offer a promising way to change experimental system 
di'sigu to computer-aided design. Important, however, is that in context of multi-body system design 
dy mimics criteria ai-e highly nonlinear and their evaluation is rather time consuming since it involves 
numerical integration of differential equations of motion. It is in the above context that we see that 
multi-objective genetic algorithms are superior to conventional multi-objective optimization algorithms. 
iN’on dominated sorting genetic algorithm -ri(NSGA-n) is capable of providing the global optima in a 



'iii.Lih* run. 

1 .3 Organisation of the present work 

1 lii‘ ])i('.s<'nt study is dividetl into five <‘liii|>tors. An outlint^ of the five <‘hai)lors is pr»?s<'iit«'d lii'low ; 

• Cliaptin- 1 : This givos a briof introduction of the present work and (iisensses tlie previtms work 
and t heir shortconiings in brief. 

• (‘hapt.ei- 2 : Tliis chapter deals w'ith the study of vehiede dynainics and snspcMision theory. De- 
sign eon('('pt and suspension principles an; also discussed. A cjualitative study of the snsjiension 
|)iiraineters a.u<l their characteristics is made. 

• t'hapter.'l : This chapter presents a study of the various algorithms that an? used in the pre.sent 
work. This also elaborates their usefulness for the pres(?nt work arul pres<‘nts a si(‘|) by step 
pi'oe('(lure for the way they are used in the present study. 

• C'lia|>t.er -1 : This ehajitcr deals with the problem formulation and its statom<mt. A dtttaiksl 
d(!Scripf.iou of the two dimensional and three dimensional models used is givtm. 

• (diapter T) : This chapter presents the I'esults of optimization of both twiwlinuinsional model an«I 
I lii <*e diiu(?n.si()ua.l model and the inference's of the rmilts are dis<;uss<!d qnalitativ<?ly. Alsu, a s<?oi)e 
for fnl.nri' work is presented and the conclusion of the whole study is given. 

1.4 Closure 

III I liis chapter w(? have given a brief introduction to the present work. A brief qualitative study of the 
lii'cvious work on the topic is also done. A comparative study of the previous work with the presciut. 
work clearly imlicates that shortcomings like the availability of a suitable multi-objective optimization 
algorithm, inadeejnate results owing to local optima, apriori knowledge of the real life problems, high 
riimpntational time etc. are overcome in the present study. This also points to suitability of multi- 
iib jectivii genetic algorithms to the multi-body dynamics problems in the engine(^ring field. 



Chapter 2 


Vehicle Dynamics and Suspension 
Theory 

2.1 Car suspension 

1 Ilf priuiary thnction of the suspension system is to isolate tlie structure, so far {is practicahli^. from 
>hofk loading and vibration diu? to iiTegularities of the roiicl .surface. Smindly. it must <lo this with- 
uiii, iiiiiiaij'iiig the stability of vehicle. The priinaiy requirement is met by the use of Hcixibh* elements 
and d.uiipi'is. while the second is fuihieved by controlling the relative motions l)(d,ween the unsprung 
masses, whei'l and axle a.ssembli<*s and the sprung masses. Tht( diameter of the tyre, size of eimtacl 
pai eli bel.weini tyre and road, the rate of the tyre acting as ii spring, and weight of whe(»l and {ixl<* 
assembly Jiffei't t he magnitude of the shock transmitted to the axle, while tin' amplitude? of the* wheed 
motion is influe'ne?etei by all thejse fae;tors plus the rate of the suspension .springs, ehunping e*tfe*e?t of the; 
shoe-k {ibsen be’rs. anel the weights eif the sprung and the unsprung masses. The; un.sprung intuss e-iui he- 
ll msely defiueel ius that beitwetem the rosid and the main smspension springs, while? tlie sprung mass is 
t hat suppeirteiel on these suspension springs, though both may also include the weights of parts eif the* 
springs juiel linkage»s. 

A systetmatic treiatment of the vehicle as a dynamic system best starts with the basic propertie-s ejf 
a ve'hicle? em its suspension system- -i.e., the motions of the body and the axles. At low freiquenenes the* 
boely whiedi is considered to be the sprung mass portion of the vehicle, moves as an integral unit em the* 
snspe-nsions. This is the rigid body motion. The axles and associated wheel hardware , which form the 
iiuspruug masses, also move as rigid bodies and consequently impose excitation fordes on the sprung 
mass. 

The dynamic behaviour of a vehicle can be characterized most meaningfully by considering the input- 
output relationships. There cma be diflerent kind of inputs. The output of interest is the vibrations of 



I ill' liudy. The rat io of ontimt and input amplittuk« reprosonts a "yam' for ll»' dyaainii- .syMciii. Tin' 
ii iiii irausiiiissihility is often used to denote the gain. Transniissibility is the nondinn'ii.sionali.scd ratio 
o| ri'spoiisi* amplitude to exeitatiou amplitude in steady .state foreed vibration. 


2.1.1 Suspension modelling 


I.iilii'ieui approaelu's for th(‘ modelling of suspension .systems are n.se<l. As an (‘xtiniple for different 
iiiodelling approaches, the vehicle's suspension chai'acteristics are considered. (l(*n<‘rally suspi'iisioii 
models can he assigiusl to three difierent categories : 


• kinematic models 


• comhiiu'd kinematic and ela.sto-kinomatic models 


• combined kinematic and ehusto-kinematic models including inertia properties. 

I load excitations arc' transmitted to the passemgers through tyres, suspension system, chasis. car body, 
and the sc'at. .\.lt hough this requires a detailed luodellingof the intersudions of each ofthe.se c'ompoiients. 
we limit our study oidy to the suspension system. Since an optimization proccxlure rec|uires iterative 
evaluation of many solutions, this simplification reduces the computational comirlexity and helps to 
.icliievc' t he optimal solution tjuieJely . If necessary, a local search with a detailed model may be performed 
in the vicinity of this solution to improve the .solution accuracy. In this study wo perform the former 
' ipi imization only. 


Sprung mass 




Unsprung mass 


Tyre stiffness K, 

Road 




Figure 2.1: Quarter car model 





A siiiiplifit'd vdiirh^ suspoiision system or a quarter ear model is shown hi the linure ‘-M. In tin 
;niiilelliii}i; oi' a suspension system, eaeli tyre is moilelled as a spring and is I'onnected to i he axle haviiej 
111 iiiisiinmg mass. The road exeitation is transmitted to tlu* axlt* throiigli tlie lyre. Kaeli axle 
I oiiiiected to the main body of the ear. which is called the sprung ma.ss. through suspension coils ami 
daiiipi'is. .-Mtlioiigh the passenger .seat is (•onm’cted to the main car luxly through llexihle connectors, 
wv sliall ignore them in this study and only concentrate on the (>ff('cl of car suspension .systems on the 
mot ion of the sprung mass. 

2.1.2 Design concept 

.\ii iiilegrai.ed design approach for dynamical .systems has to supiiort all steps from problem formulation 
In prolilein .solution by (jptimixation. Firstly the technical system to be optimized are transfornKsi to a 
mat lieiiiatical moflel. The mnlti-body system approatth is widely used in vehicle dynamics. Then, design 
goals are dc'fiui'd which is often difficult, since technical requirements and human wishes tire sometimes 
ilillii ult to be formulated fus matiunnatical functions. Besides the complexity of tin* mtidels this is maybe 
niie ( if I he reasons why even integrated design methods cannot stibstitute a design engineer, but supiiort 
b\- sol'tware systems helps to make design better. In order to improve technical systems, design changes 
have to lie nuide. Therefore, parameters of the model are classified either as design vari.ihles whose 
values can be varit'd for optimization process or as system constants whase values are fixed during 
opt iiiiizalioii. In our sttuly we consider the sprung mass, unsprung m(us.seH. wheel to wheel distances 
as the system constants and the stiffness and tlie damping parameters tis the design varitihles. After 
all these preparations the optimization part of the design proc(*ss may be started. In ease of .severftl 
uli jectives. a multi-objective apprtjach is applied. 

2.2 Suspension theory 

When ti rigid body, such as a car body, is mounted on springs at both ends, and the front .springs are 
siibji'cted to vertical forces through the front wheels, oscillations are produced in both front and rear 
springs. A similar action occuns when the vertical forces are applied to the rear wheels . One of the 
ea,rli<!st treatments of this reaction between front and rear suspensions was presented by Professor .1..]. 
(luest [12]. 

For .simplification, lateral forces are not considered. Moreover it is also assumed that both front 
\vh(!(!ls are subjected to identical vertical forces simultaneously. The same conditions of identical forces 
also apply to the reax sprinp. The mathematical model is, therefore identical to two pogosticks con- 
iicvited by a rigid bar as shown in the figure 2,2. Guest’s model is shown in the figure 2.3. If a vertical 
force is applied at the front end of the beam, the line Xy which represents the center line of the body, 




Figure 2.2: Simple coupled .suspensions 



Figure? 2,3: Guest’s model : Elastically conjugate points 

will !idoi)t an inclined po.sition X'Y' , rotating about point B. If a vertical force is applienl at the* re'ar. 

I lieu rotation will be about point A. Any pair of points that possesses this relationship can l)e shown to 
sat isly the e<piation pij = ab. The dimensions a and b are established by the relative rates of the IVoiil 
and rear springs ki and since kia = hib. The spring center C is, therefore, the point of balance 
and a vc'rtical force applied at C will lift the beam with no tendency to tilt. Points A and D an? called 
{ laHtimlly conjugate points. 

'riie next stage of the Guest construction [12] is the dynamically conjugate points of the sprung mass. • 

The moment of inertia of the sprung mass about the centroid is MK^. The sprung mass M for which 

I,, = MK'^ can be represented by discrete masses mi and m 2 at distances r and s from the centroid G. 

'I'o satisfy the equation, ra = the following conditions must be met: 

mi + m2 = Af, 
mir = m2S, 

mir^ + m2S^ = MK^. 

As in the case of the elastically conjugate points, the quantity r is chosen arbitrarily and the 
corresponding value of a is found frorp th« rdaiioasfaip rs = if®. Just like the elastically conjugate 



!)iiiiii.s then? an' iufinito iiunilxT dynamically conjugate points. 

riii' final stage is to comhine the two syst(nu.s. This will give ns the doiihlr nmjuiiuh [niiiii.s iliai 
w ill satisfy bot h conditions. There is only one such pair of points. 

To i-ahadate double <a)njng<i,te points, the following information is us<'d. 

pq = = nb. 

r.s = A'“. 
p = r + ic. 
q = .s — a:, 

wliere is the distau(t<‘ between the center of gravity 'G' and the spring ceiittn- . a ami h an> the 
(listan(a\s of the spring (-enter from tlui front and rear wlieels . L is the distances betwwni the front and 
liu* rear axh's and k\ ami k-> are the front and rear spring rates nsspectively. In figure 2.4, G is the 



Figure 2.4: Double conjugate points : A simple mathematical model 

centre of gravity, C is the spring centre, and H is the conjugate centre. 

K’^ is a (juantity that is measured experimentally. Usually it is taJken as 




where li and I -2 are the distances of the center of gravity from the front and rear axles. 
Prom the above equations we get, 

K'^-a^'-eA J abf + 

r SS — r 4- ' ' ■ - — — ^ 

2 ® 2 ® 



,\ few Special <-a,s(W av(' discussed balow : 


• Mi)St ni(td('rn vc^hiclcs with substantial front and rear ov<»rhang exhibit <lynanti(' index clusc in 
unity. Tlial is : 

DI = K'^/lih = 1 

WIk'u till' ('(piality holds, tlu* front and rear snsponsknis are lo(‘af.ed at conjugate c(“nt('rs of per- 
t ussiou (an input at one suspension causes no reaction at the other). In tiiis avtn' ti«> o.scillutioii 
centi'rs an' located at the front and rear axles. This is a desirabh* condition for good rid<! as tliere 
is no interaction betwei'ii the front and the rear suspensions. Some special castw are di.scu.ssed now 


• Spring c('nt(!r at the CG This condition corresponds to the tuise of couphitd pitch and Imuuce 
motions. The ]}itch juid bounce oscillatiojis are totally independent. Poor ride results becausi' the 
motions can be very irregular. 


Dynamic index K greaten- than unity This occurs when there is substantial overhang at the front 
a.iul/or n'iir end of the car. Th(* bounce center is in front of the CG (beyond th<‘ front tixle) and 
t he pitch <'(!nter is betw(H*n the CG and the rear axle. The natural freciueucy and tltit. ride will 
still result if t,h<^ spring center is located far tmoiigh behind the CG (front ride rate h^ss than rear 
ride rate). 


• I’mumph'd motion and dynamic index equal to one- This cumlition results in iHpial bouiu-i’ and 
pit ch fri'cpiencies. Tlu! ride is inferior because there is essentially no pattern to the road- geimrated 
motion; it is quite unpredictable. 

2.2.1 Suspension elements 

The snspc'nsion elements present in the model are the springs and the dampers. The suspension spring 
is in series with a relatively stiff tyre spring, so the suspension spring predominatiis in establishing the? 
rifie rate and, hence, the natural frequency of the system in the bounce (vertical) mode. Therefore, in 
ihc ])resent study we have considered tyre stiffness as a system constant while the suspension spring 
si itfuoss is a design variable. Since road acceleration inputs increase in amplitude at higher frequencies, 
i.lic best susptjiision is achieved by keeping the natural frequency as low as possible. For a vehicle with 
a given weight, it is therefore desirable to use the lowest possible spring rate to minimize the natural 
frequency. A softer spring leads to lower acceleration but the spring stiffnesses are usually kept above 
(tertain vadues to make the ride handling easier. 



D.iuipins av<i list'd to mako tht* viliratious die down. In alisoiifti of daniphif!; a vibrafiou oiifc cxfiiod 
w.mld t'xist iiidt'Hnitoly. In pmcticc damping exists in three forms : friction, viscous and lliat due to flic 
pri'st'iicc of air. 'I’lic suspension of most vehicl(*s exists to keep to a minimum the disturhaiict's of the 
sin iitig body due to irrt'gnlarit ies in the tt'rrain over which tlit' vehicle runs. Friction, in the su.s|»ension 
liiilcige. me.ans that no movement of the .suspension itself occurs IVir irregularities below thost' .at whicli 
I iie dt'Hectitm of tht* tyre imjdies an increase in force etinal to static friction force. 'I'liis constitutes only 
;i small component of the total damping. Damping due to air is also small compared with that necessary 
lu pri'vent unduly Large movements of the sprung weight and can equally he regarded as insignificant. 
The main Imlk of ilanqiing is therefore likely to be viscous damping from the hydraulii; dampers whicli 
are now iniiver.sally n.s«'d. It is referred to tis viscous damping and taken as bi'ing projiortional to tlie 
Viatical velocity lietween unsprung and sprtmg masses becaasi' this is mathematically true and is rea- 
soii.ihly clo.s(' to till' true state of affairs. 

The amount, of dam{)ing required is, as in so many cases in engineering, a compromise ; in this 
case bi'twi'en that needed to jtrevout undue persistence of a vifiration at natural frequency excited by 
a single disturlianci'. which also prevents the build-up of exce.s.sive amplitude of a forced vibration due 
lo a series of impulses, and the uncomfortable fact that the higher the dfunping fon-e the greati'r the 
disturbance fed into the .sprung nuuss by any given road irregularity. The damping force is resiionsible 
for an additional distnrliance on the sprung weight over and above that due to tin* <lefIe(;tion of tin* 
suspension spring. 

Of the two different frequencies : one Is that of the sprung mass on the suspension spring system. 
;iu<l t.lie second is that of the unsprung mass- the wheel and axle assembly- on the tyre. Obviously th<‘ 
latter is affected by the suspension spring rate, but only marginally. The former will be experienced as a 
relatively low fnuiuency- perhaps about 1 to 2.5 Hz-bouncing of the carriage unit, while the latter is that 
o| wlu'i'l Iiop. at a higher frequency-generally 10 to 15 Hz--and is generated almost totally independent 
■ if tlie motions of the carriage unit. For minimizing the amplitudes of wheel hop not only re.sonant 
but; also isolated hops - the unsprung weight must be kept as small as possible. Resonances of either 
t he sprung or unsprung masses can affect adversely, and indeed to a dangerotos extent, the handling 
characteristics of the vehicle. Obviously, therefore, it is important to maintain the dampers, or shock 
. ibsorbcrs, in good working condition. 

2.2.2 Suspension principles 

Pure bounce will occur if the front and, rear sprung massif are equal [14], the front %nd rear springs 
have identical frequencies and identical rates and arc in phase. A car can be designed to have a 50/50 



ui'i^lii (listrihntiou and identical springs all annind. but tho pluusing of ripples in the mad surfarc an* 
iii'vimd onr control. It is iiu'vitabk* that a dt'gnto of ijitching must oftc'u occur. ’J’hc inertia of the 
spiling mass will resist this pitching motion and the dispo.sition of the major coinponcnt.s. i.c. the 
t'ligiiic. i.ransniission. fuel tank, passengers, etc., that make iij) the sprung mass <‘ontribute to the re- 
-isiaiici* ('xerti'd by the sprung in opposition to th<? pitching moment. If the major masses of the body 
1 ompoiieiils tend to be n(‘ar tin* <‘ent.roid. the sprung nuuss offers less resistance* to pitching than a body 
o| ilie saiiK* total mass t hat has the more nuissive components situated at a greatin* distanca* frotn the 
ci'iiiroiil. Expjt's.sed tt'clmically. we say that the second case lists a Isirgcn* polsir nuiment of iiu'rtia. In 
ihe present, study we are intitrostetl in two polai’ moments of inertia, i.e.. pitching moment of im^rtisi 
, 01(1 rolling moment of inertia. .Although utany modern cai's hsive the centroid well forward, the value of 
; l\ - is often close to unity. This is more by chance than by deliberate design. K is the radius tif 
gvrat ion and l \ and / j are the distanetts of the center of gi’avity from the front and rear axlits reK])ectively. 


A gi'iieral rtile tliat is observed in a vehicle suspension design is that the rear natural fretiuency 
,<iionld be greater than the front natural frequency by at least 10 percent. This is rationalized by the 
( ibsei val ion l lial vehi(;Ie bounce is less annoying as a ride motion than pitch. Since (?xcitatiou inputs 
IVoiii t he roa,<l t.o a, car affect the front wheels first, the higher rear to front ratio of frequencies will 
lend to indmu* liomux*. Also the locations of the motion centers are depeiuhnit on tlie relativi* values of 
I lie natural fr(’<inencies of the front and rear suspensions. With equal frequencies it is found that mu* 
(•('liter is at tin* Cl.CJ. location and tlm other is at infinity. Equal frequencies correspond to decoupled 
vcrtic.il and j>it.cli modes, and "pure" bounce and pitch motions result. With a high<n’ front, fretpieney 
I lie motion is coupled with the bounce center ahead of the front axle and the pitch center toward the 
r(;ar axle. A lower front fretpiency puts the bounce center behind the rear axle and the pitch center 
forward near the front axle. This latter case was recognized by Maurice Olley in the 1930’s [12] as the 
best for atdiieving good ride. 


Since the suspension spring in series with a relatively stiff tyre spring, suspension spring predominates 
ill estahlisiiing the ride rate and, hence, the natural frequency in the bounce (vertical) mode. Since road 
acceleration inputs increase in amplitude at higher frequencies, the best solution is achieved by keeping 
t he natural frequency as low as possible. For vehicle with a given weight, it is therefore desirable to use 
I he lowest practical suspension spring rate to minimize the natural frequency. However, the practical 
1 ill ill, s of stroke that can be accommodated within a given vehicle size and suspension envelope constrain 
t he natural frequency for most cars to a minimum in the 1 to 1.5 Hz range. Performance cars on which 
ride is sacrificed for the handling benefits of a stiff suspension will have natural frequencies upto 2 or 



.! 'i 11/. Ill tlic ]m’S(!Ht stiuiy. all flic I'miiuniiut's i.o. frout. rear, pitrliiiifs; aro runslraiiiml in a ramj;t' ni 
II. s In 2."> Hz. This allows for a hot, tor handling of a car hosidos inaintaiuiug a good ride fniufori li'vd. 


2.3 Rigid body bounce/pitch motions 


Thr .siiii|)l(‘ iiii'chaiiics of tho tiiiartor-i’ar model do not fully reiireseiit thi' rigid-body mol inns lhal may 
. Hciir nil a motor vohiclo. Because of the longitudinal distance between th<' axles, it is a mult i-ini>iil 
systi'iii that responds to the pitch motions as well as vertical bounce. This is the reason for tamsitlering 
I he two dimensional model instead of a quarter-car model in the present study. Dejiending on the 
mad and speed conditions, one or the other of the motions may be largely absent, or they may not 
necessarily be ol»serve<l at the point on the vehicle where the vibration measurements are made. The 
jjiich motions an' important because they are genei-ally considered objectionable and are the primary 
source of lougitiidinal vibrations at locations above the center of gravity. Understanding the pitch and 
bounce motions is essential because it is their combination that determines the v(?rtical anrl longitudinal 
vibrations at any point on the vehicle. 

2.3.1 Interaction of front and rear suspensions to single applied disturbances 


.\ car may meet a single' bump or hollow affecting one track only, or it may meet a traiisver.se ridge or 
hollow afieeting both tracks. If we consider a transverse ridge, as a result of passing over it the front 
end of the car is raised. Then, after an interval depending on the wheelbase and the spi'ed. th<' liack 
end of the car is raised. 


TIk' lirst effect of lifting of the car is to excite some pitch movement. The subseipient lifting of the 
liaek end of tlu' car may mitigate or aggravate that pitch movement, depending on the time interval 
and the pitch frequency because the movement of the sprung weight on its suspension is not critically 
(lauiped. any movement of the suspension at each end of the car will persist with diminishing amplitude 
i'lir a iiiuuber of oscillations. Depending on the relative frequencies of the front an<l real’ suspensions 
and t lic amount of damping, the initial pitching motion may be reduced or increased by the passage of 
1 iiiu'. A few facts about the above case are listed below : 

• With the lo'wer frequency rear suspension there is a strongly persisting pitch movement. 

• With the higher frequency rear suspension the initial pitch movement is very quickly reduced. 

• With the higher frequency rear suspension the initial pitch movement is never so severe as with 
the lower frequency. 



• rii(' initial ])it.ch inovoniont is loss sovoro at the higher than at the lower speed, la-eaiise ilir liiiif 
iiit(*rviil is l('ss. This tinu* interval is at both spetnls less than half a cycle of the front siispensinii, 

2.3.2 Effect of regularly repeated disturbances 

Many rotid surfaces have develoj)e<l fairly rt'gnlarly oc<airing bnnips and lutllows. For a car traversing 
siieli a ro.wl the wheelltase filter effect has to b(‘ con.sidered. In its siini)lest f'onn. for condif itm wliere the 
whicle wluH'lbiise etpials half the road disttirhance wavelength, pitching is exc.itttd if the conibinjil ion of 
\'eliie|(> speed and disturbance wavelength produces a disturbance freqtiency equal or trlost* to tlii^ vehicli'‘s 
pitcli friHiuency. The ('xtent to which the pitching amplitude builds up dttpends on the suspension 
.l.ini|>ing fiictor. Somt! build ui> will also occtir if the wheelbase is 3/2. 5/2 ett\ times (he disturbance 
\v.i\el<(iigl li and t he vehittle .spetid productts impulses at or clostt to tint vtthicle pitch fretUHtiicy, This 
elleet eousiderably affects the vehicle's vertical acceleration spectra at the cviitKr of gravity find rc'.-ir 
wheel positions, find the pitch acceleration spectra. 

2.3.3 Relative pitch and bounce frequencies k^/l\l -2 ratio 

Assinning that tin* frtsim'ncy of oscillation is unaffected by pitching fis opjfosed to bonncing inovisneiii 
III till' c;ir. this is generally nearly true. To be true the front and ntar mid imisses would hfivt* to lie 
l•l^e(•tiv^'ly coiKUJutrated in front and rear axle vertical planes across the cai-. If this <‘oudition is met 
I hen the moment of inertia I of the car about a transverse axis through the centre of gravity as shown 
in the ligtin! 2.5. is : 



Figure 2.5: Mass diagram, = 1 

/ = W}1\ + Wrll 

because a and b are the distances a£ the centre of gravity from the front and rear wheels 

Wfh * Wrh 



Sul)stirul iu}j: ill tlu' sc'cnnd term of tho monunit of inertia ('(luatioii. wt* gi'l : 

I = Hy, (/, + i.) 


Xiiw. 11'/- = Wltj L . HO 


/ = U7,/a(/i +h)IL = W'hh 


riif usual expression for the ini'rtm / of a body is WK^ where K is the radius of fi;yratinii. It follows 
liien'fore that tin? ratio as an indication of how far the car departs from the characteristic of 

i'(|nal jiitcli and bonni^e fretiuencies. 


riie fact that K-/I.\l .2 is {■•enerally Usis than 1 has two effects : 

• TIk' iiitch freciiiency is higher and therefore for given amplitudes the accelerations will he greater. 

• rh'caiise the resistance to pitching is less a given exciting force will produce a greater aniplitnde of 
jiitirh. Looking at a vehicle suspension with tyre acting as a spring, the less the resistance of the 
sprung mass to the pitching motion resulting fi'om a given road ridge as disturbance the less the 
tyre deflection and therefore greater tlio deflation of the vehicle suspension system ami therefore 
t he greater the initial pitch amplitude of the sprung mass on its suspension. 

ill ideal case K“ /I \l) = 1 is used in the present study and the centroid is nearer to tiie front tixle. so 
iliis (;onihiuatimi lemis to the faster dying down of pitching motion as it is the most irritating motion. 

2.4 Non-linear suspension characteristics 

.\on-linearity is nbicpiitous in nature [15]. Linearity is an approximation to reality. In shock and vibra- 
lioii syst.emH, isolators such as air springs, elastomeric dampers, and wire-rope isolators are inherently 
iionlimvir. Assumption of Hooke’s law for springs and linear viscous damping for dampers is done just 
lor mfil.heimit.ical simplicity. Sometimes the amplitudes of steady state vibration are small enough l.o 
Justify the assumption of linearity. However, the transient displacements may often be sufficiimtly hu-ge 
wlieu the non linearity in springs and dampers cannot be ignored. 

It has been established that the effect of nonlinear damping is more predominant than the effect of 
nonlinear stiffness so far as the performance of shock isolators for base excitation is concerned [6]. In 
I he jnesent study a cubic non-linear damper is considered in addition to the linear viscous damper, in 
( he case of shock excitation of the suspension system by a rounded step function. However, the stiffness 
characteristics are considered to be linear throughout. The damping force expression with non-linear 



is .u;iv(,'n as : 


F = c.fd + Cf„(B 


2.5 Closure 


111 ili(' pivsi'iit ch.-iptor wo (Usaissod th<^ car .susj)eij.sioii iiiodollinK and il.s dosigii . lu im.lor, stand ilio 
pivsciit woi-k i)n)p('rly a iwiofdiscnssion ou t,ho susponsion thwiry is uiadtt. Thi' snsponsicm |)i in< ipl(>s ate 
als.) incnlionc'd as Mioir kuowlodge is a unist for the problem forimUtition. The suspension pjinmi(‘t<‘is 
aiv also studi(>d in bru'f and tlnnr effwtt on the ride characteristics is also inentioiK'd. A (imilifatiw 
..I ndy of li'Iai iunship between pitch and the bomice motions is made and the efietd, of dynami<- index 
. ai t hcse mol ions is also discussed. At the end of the chapter, non-linearity in the damptn-s is dhscnssed 



Chapter 3 

Algorithms Used 


;i.l Introduction to genetic algorithms 

( Algt)rilhni.s (GA) are finding extensive use now-a-dayvS in the field of optimization relat(xl prob- 
lems. They form a V(tv robust tool for the optimization pai't of computer aided design process. Gtnietic 
.ileoritlmis ai-e more suited for optimization when used alongwith a numerical integration technitpie. 
Ill i li(' pnwc'iit study we have used different kinds of genetic algorithms, i.e. real codcsi gelK^ti(• algo- 
riihiii and the multi-ol)j(H:tive genetic algorithm.The genetic algorithms otter many advantage's over the 
cniiveiil ioiial methods of optimization. A few of the adviintages are : 

• I'nlike convemtional technique genetic {Ugoritluns ai'e robust, i.e. they can be useei t'ifie-iently for 
a wide range! of j)re)blems. 

• GAs are' not senisitive te) initial solution chosem. 

• G.\s elo ne)t reepiire! eleirivative or any sue:h information. 

• GAs are* e;apable of finding global optimum in just one run. 

3.2 Standard genetic algorithm 

( ie'iK'tic Algorithm! (7 j 4) is the optimization technique, which mimics the nature in order te) finel the* 
global oi)tima of a given problem.This is inspired from the concept of “Survival of Fittest”. It meieuis 
1 hat, only those individuals, which are better than others in one or other aspect survive. Unlikej clas- 
sical optimization methods GA deal with a population of individuals. In Binary GA individuals are 
rcpiesented by the binary strings called chromosomes. These chromosomes axe the string of O’s and I’s 
callt'd allele. -In Real Coded GA, chromosomes comprise of string of real numbers. Salient features of a 
standard simple GA are as following. 


1. Selection. 



‘J. {,’ross-uv('r. 


.*». Mutation. 

:}.2.1 Selection 

i'cir a si'cii poimlation <tf iiulividuals only thoso imlivulimls an* solcctoti wliicli an* Ix-tfor in li(iM's.s than 
I ii liiTs. This way. G.V follows tho Darwinian law of survival of t h»' fittest an<l hy this sclcrt ion our scarrh 
is i^iiidod towards th(' optima. More importantly, this mothod does not rcjoct the solution which is not 
i licsi. I'his .solution may have tho information, useful in finding the global l>ost sohition. Most used 
iik'iIkmIs for scks'tion are 

1. ToiirnanK’nt Solo<-tion. 

2 . Stochast ic Universal Sampling. 

Ronh^ttc Wland Selection. 

;{.2.2 Cross-over 

111 cross-over two individuals mate and prodtico children. This is a very important operator which crctitc 
I he new solutions, by using tho [woporties of the j)ai*euts. This create the solutions in tho region «»thor 
I halt th(’ t'xisling region. This explores the search space to find the best region. It is supposo<l that the 
iiulividuals which are surviving, have some good string, which let them alive. So good parents combine 
lo pi oduc(> botti'r children. Two types of crossovers ai'C mostly ustid: 

1 . liimiry C’ros.sovor. 

Deals with chromosomes with binary numbers. 

2. Real Cro.ssover. 

D('als with chromosomes with real numbers, 

:i.2.3 Mutation 

riiis mimics the Hud<len changes occurring in the nature, that creates the individuals, different from the 
individuals of its own category. This operator is particularly helpful when the search has stuck at some 
local optima or search is going slowly. This can also create solutions which are not good. They will be 
taken care by the selection and other operators. 

Tho algorithm of simple GA is shown in figure 3.1. 




Figure 3.1: Siiopte genetic algorithm 








:i.3 Real coded GA with SBX operator 


|}iii;iiy (-(xled GAs discreti/o tlu' st’umrh space by using the coding of real variables in Itinarv siiiiig>. 
However, the search space is real and ctnitiniions. The coding of ri'al variabh's in finite-length si rings 
r.nises a niiniber of probhans lik(’: 

• Inability to achi<’V(' arbitrary pnH:ision in obtained solution. 

• Fixed mapping of problem variables. 

• Inhen'iit Hamming Cliff problem Jissoiuatod with binary eroding. 

• l^riaressing of Holland’s .schemata in continuous search space. 

lb ovi'i-come all such difficulties the real coded GA (RGA) are used. In the present work, for single 
iihjeetive optimization an RGA using Simulated Binary Cross-over or SBX [2] has las'ii usisl. This 
l !C!.\ (unploys a real-coded cross-over whose search power is similar to that of the single point (rross-ovin- 
used in Ijiiiary-cisied GAs. In order to achieve that we define the search power of a cross-over operator 
ill tenns of t lu' jirobability distribution of any arbitrary child .string to be created from any twti given 
ji.neiit strings. The .search power of the single-point cross-over operator is first calculated. Later, the 
>imtil;i1ed binary cross-over (SBX) is developed having a search power simihu' to that of the single-])oiul 
I Toss-ovt'r. The diff<w<;nce in the implementations of the real coded GAs with SBX and hinary-cotled 
( lA.s with singhvpoiut cross-over is that in tint former method the coding of variables is elimitiated and 
I child St, ring is (treiitt'd from a probability distribution that depends on tim location of th<; parent strings. 

Ill tlie <;ase of simulated binai'y cross-over (SBX) a parameter called '’’spread factor", denott'd by li 
is detiiied. Based on this parameter the cross-over operator can be of three types. Operators having 
i' < 1 nre called contracting operators and those having P > \ are called expanding operators. The 
uperators having /f = 1 are termed as stationary operators. The random variate ^ is generated using 
I lie polyiifimial distribution as below, 

p()8) = 0.5(?7 -t- for yS < 1 

p(yS) = 0.5(ij + 1)-^, for ^ > 1 
77 - 1-2 

Th(^ value of ’’spread factor” is calculated by equating the area given by a randomly generated random 
number between 0 and 1 to above equations. The probability distribution curve is shown in figure 3.2. 
iiere t/ is the distribution index of the probability distributfon. After calculating from above equations 
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J2 



Fii»;un' .‘{.2: 


SBX ojxM-ator. probability distribution for thr gt‘noratx»n of ohildron 


! Ik* cliildrtMi points aro )i!;oii<>ratod as follows : 


y} = 0.5[{H- A)a:| + (1 - Mxf] 
lif = 0.5[(1 - 0M + (1 + 


I i I'll' is (U^iiiu'd cis 


■ii; 


|. The sprcful fac*t.or is rancloinly generated from a polynomial ilislri- 


I (III inn (lelined above for cadi variable in case of “variable by variable” crossover. 


The ji is generated l)y transformation of above polynomial distribution using nniform variate fis 

Ii (111 (\VS. 

Ii(v.) = {2u)'ii~>, if «(0,l)<0.r) 

f3{u) = [2(1 - if u(0, 1) > 0.5 

II. is a uniformly generated random number between 0 and 1. 


In ('iisi’ of multi-variable problems, the crossover is done “variable by vai'iable” i.e. for a variable the 
crossover will tiike plac(J with a probability 0,5. If the crossover is to be done the value of ft is cal(.mlated 
using above equations. 

The chai-ac.teristic of the SBX operator is that it generates offsprings depending on the position of 
I lie [)an'nts. If the parents are far, the children solutions are created fax away and if the parents are 
close the children are generated close. Also a small value of the distribution index creates points far 
away while a larger value gives closer solutions. 

Simulated binary crossover for problems having no bounds on the variables has been described above, 
ff the problem parameters have bounds, SBX can be heed alter some modification. In this ca^ the area 



ul i 111 ’ in iihahility distiihutiou which lies onlsiih' the Inmmis is iuhUii fu t ii.il is ciicliiscii Itv i he iMuiiKi- 
Mailu’iiuilic.-illy it is done as follows. 

L_ 1 

= (Hft-) . it «< — 

o 

. , 1 ';H, , 

f^i/ = n )• otlimvisc 

2 - H,a 

wlii’i'c. o = 2 - and ii is calrulatcd its 

2 

,:< = 1 + - .?:/) . - To)] 

T2 - .P| 

I Icri'. XI and x„ are the lower and upper hounds on the [)arameter ,v. 

3.4 Multi-objective genetic algorithm 

ill a classical oi)tiuii/iatiou approach only a single objective function is minimized leading to a unique 
solut ion in geinn-al. Pra(;tical design problems do not look like this. Dynamif' systems have to he oi>- 
liniized with I'esjxnrt to several conflicting specifications and several different designs an^ acceptahh' as 
opi iinal with i-espiH-t to th<* same set of specifications. Multiple optimal solutiotjs exist because no one 
solni ion can he optimal for multiple cxmflicting obje<;tives. Once such multiple solutions an* found, tisn- 
ally. a higher-levfd de<usion-making strategies is adopted to choose one solution from tin? set of obtained 
i’aii'fo-uptiinal solution.s. 

The principles of nudti-criterion optimization ai’e different from that in a single-objective optimiza- 
! inn. 'file main goal in a single- objective optimization is to find the global optinnil solution. However, 
ill a mnlti-cri1;eria. optimization there are more than one objective function, each of which may have a 
ililferent individual optimal solution. If there is sufficient difference in the optimal .solutions corr<ispon<l- 
iiig to different objectives, the objective functions axe often known as conflicting to each other. In the 
present study the two conflicting objectives considered are the minimization of vertical displacement 
and minimization of vertical acceleration . The classical ways of tackling multi-objective ()i)timization 
inobletns is straight forward (16, 17, 18, 23, 25]: Convert multiple objectives into one objective. There 
I’.Kists a number of conversion method- weighted sum approach, c-perturbation method, Tchybeshev 
method, inin-max method, goal programming method and others. Since multiple objectives are con- 
vei tiid into one objective, the resulting solution to the single-objective optimization problem is usually 
subjective to the parameter settings chosen by the user. Moreover, since usually a classical optimization 
method is used, only one solution cam be found in one simulation run. Thus, in order to find multiple 
I’areto-optimal solutions, the chosen optimisation a%6rithm must have to be used a number of times. 



l intlK'nium'. tlu' classical niot.hods have been found to be sensitive fo the convexity and l•olllinuil> >i| 
I lie I'.uctn- optimal region. In this study a multi-objective geneti<; algorithm is used in conjimciion 
vviili a mmierutal integration tia^hnujue. 

.\hilti-objectiv(' evolutionary algorithms with such conflicting objectiv<> functions give rise to a s.-) 
.il opliiual solutions, insteaxl of one optimal solution. The retvson for th(^ optimality of many solutions 
ir. ih at no oiu' can he considered to be better than any other with respect to all objec^tive functions, 
TIu'se optimal stilutiuns have a special name Pareto-optimal solutions. 

The concept of oi)tima.lity in multi-criterion optimization deals with a number (or a set ) of solut ions, 
instead of um\ sohition. Bfused on this, we first define conditions for a solution to bwtome doininal<>d 
u iili respect to another solution and then present conditions for a set of solutions to InTomti a Pareto- 
< \\)\ iiua.1 


I 1 I (irr.fo^ opt WKil Front is tha sot of solutions which (ivo ayuolly pood as cojupawd to any 
iithcr .•^(dul.hv of the. Hvt i.c. any aolution on the front can not be selected when compared with other. 
irilhoHl hvinii biased towards other. Mathematically it CMn he given as : '‘A decision vector x € A'/- is 
said to he non-dominated regarded to a set A C Xf if 

^ae A: a>- x 

If it is clear within the context which set A is meant, it is .simply left out. Moreover, x is .said to he 
!' an to- optimal if x is non dominated regarding Xf.” 

i' igiire :{.4 gives a bettor illustration of the concept of Pareto- optimality through a example problem ijf 
n ad<>-ofl' l)etw<^on two (conflicting objectives namely, cost and the accident rate, both to be minimiztsl. 
I'lie [joint A represents the solution with a very small cost but highly acccident prone. Another solution 

II l•('pr(^seuts the .solution very safe but very costly. Here both the solutions are good with respcict to 
t he. one objtKctive, but worse in other objective. There exist a lot of solutions (like solution D) which 
also Ixcloug to the Pareto-optimal set and have equal importance to a unbiased decision maker. All 
i luisc solutions (lying on the thick solid line) are known as the members solutions (lying on the thick 
solid liu(i) are known as the members of Pareto-optimal Front and this solid line is called Pareto-optimal 
Front. 





Accident rate 

A 



Cost 


-Pareto-Optimal Front 


Fillin' 3.3: Tht: concept of Piireto-optiinal solutions is illustniti'd. 


3.5 Dominance 


This is ;i v('ry iiiiportant concept so far as the discussion on multi-objwTive problems is <'onc(Tned. This 
liu iiis I he back bone of the ranking strategy. For constrained and unconstrained dominance defiuitu)n 
v aries. Hot h of tluv definitions are given as following. 

3.5.1 Dominance for unconstrained problems 

Dfifiuition 2 If their are. N ohjeeJ-ive functions and two individuals, “a” and ’‘h'' are compared then 

a >- ft (a dominates b) if f (a) > f{h). 

f{a) > /(ft) means that Vsi € N /,(a) > fi{h) 
or 

a > ft fa weakly dominates b) iff (a) > f{b). 

/(«) ^ f{h) means that 
^ieN fi{a)>fiib) 
and Ai e N fi{a) > fi{b) 
or 

a ~ ft (a is indifferent to b) if f {a) ^ f{b) 

A /(ft) 1 f{a) 


III simple words it can be said, that “a dominates b" if and only if for all function values, the function 
value of a is either better or at least equal to that of function value of 6, and at least one function value 
I if a is better than ft. If all function values of a are equal to those of ft or if in some function values a is 
1 letter and for some function values 6 is better, a and 6 are said to be indifferent to each other. 

I'liis can be understood with the help of the Figure 3.4. It is observed that the solutions which are not 
(111 the Pareto-optimal front (like solution C) also exists. If it is compared with solution A, according 
to definition 2 both are indiffereiit. But if it is compared tyith another solution D, it is found to be 
dominated by solution D. Since for -a, solution C, there exists a sedution D in the seardh space, wiad:i 
dominates it, this solution cann’t be in Pareto-optimaJ. Fiont. 



.5.2 Dominance for constrained problems 

1 iho limit (;aso.s where the coiLstraints are present, tieiiiiitiun of "(hiuiiiiaiici" is rnleliiieii in 

[•i;uuni.u)i]a,t(' the eou.straints. 

)efiiiition 3 The. definition of dominance whan nioiv. then one aonHtraints arc prcucjit i.s i/irni ai-nnilinti 

I) /./;(' Jollowinij condiiionH. 

.’oiidition 1: If in dividual. A is feasible and individual B is not then .4 dominates ]J. 

illondition 2: If individual B is feasible and individual A is not than B dominates .1. 

Condition 3: If both of the. individuals are infetisible then the. individual with lesser violation of con- 
st raint(s} dominates the other. 

rtoiidition 4: If both, of the. individuals air. feasible, or equally infeasible, dominance is decided (i< <-orilinii 
to the. definition 2 

There are priimu-ily two goals that a nmlti-criterion optimization algorithm must, t.rv to achieve: 

• Gnidi' (he search towards the global Pareto-optimal region, and 

• .Midutaiu ])oj)ulatiou diversity in the Pai’eto-optimal front. 

3.5.3 Constraint handling 

I'll!' steps for constraints handling are given as following: 

1. Xuniialise tlu'jerk and the frequency constrfunts. 

2. I'Mud tlie .square sum of the constraint violation (if any). If none of the constraints is violated, 
corresponding violation is set to zero. The function is considered to be equal to sum of it’s value 
and the square sum of the constraints mulitplied by a penalty parameter. This is termed as 
penalty function method. However, this leads to distortion of function, so it is not preferred. In 
till' present work we are using tournament selection scheme for the constraint handling purpose. 

3. Tournament selection is done to handle the constraints as suggested by Deh, [20, 19]. Tournament 
Size taken is two. Winner of the tournament between two solutions is decided by following the 
constraint domination conditions. 

• Condition 1: If any of the individual is feasible and other is not, select the feasible individual. 

• Condition 8: If both individual are infejisible select the individual with least error or con- 
straint violation. 



CcmdUion :i: If both individuals aro feasibk* or both an' iiih'asiblo with .same «t»ii>iraiiii 
violation, seleot tho individual with better fitness. 


.6 Srinivas and Deb’s non-dominated sorting genetic algorithm (NSG A 


Ik' idea behind the NSG A [21] is that a ranking stdeetion method is !is(h 1 tti emphaKiz<* currciii imii- 
oniiiiatcHl points and sharing funetion method is used to maintain diversity in the population. .\'S<!.\ 
iifers I'roiii a. simple gc'uetic algorithm only in the way the seleetion operator is us<^d. The eross-over and 
nutation opc’iators itniniin as usual. Before the selection is pertbrmed. two i)ro(;edures are perforun-d 
I'liaHv: 


1. Classifying a population according to non domination. The fii-st set of non-domiuat<«l s«ilutioiis 
> .isstuiKsl to (amstitute tin; first non-dominated front in the population. In order to find tlu* solutions 
ii'loiigiiig ttt tlu' second level of non domination, the solutions of the first levttl of non-domiiiiitu»n cjiii 
je hnuporarily disntgarded. This procedure is continued till all the solutions are c.hLssified into a level 
if Hori-dominatiou (rank). 

2. Filiiess assignment using sharing function: Initially the fitness is {tssigned to emb level according 
hi iis rank. .\n individual in a higher level gets lower fitness. This is done in order to maintain a 
selection pi-cssun' for cluxjsing solutions from lower levels of non-domination. Tims, the search is g\iid(‘d 
Inwards th(' Pa.reto-<)ptimal region. To obtain diversity among the current nondominat«'d solutions . 
lonely solutions arc emphasized using a sharing function strategy. So the fitness assignment is done in 
I wo stages. Firstly, all the solutions in the first non-dominated front arc assigned a fitness (dummy 
lit ness) e(|UHl to the population size. Based on the sharing strategy, if a solution has many neighbouring 
solut ions in the same front, its dummy fitness is reduced by a factor and a shared fitness is computed, 
riiis factor (usually known as niche count) depends on the number and proximity of ueighbouriug 
solut ions. Thereafter, individuals in the second non-domination level are all assigned a dummy fitness 
f({U!i,l to a number smaller than the smallest shared fitness of the previous front. Then, shared fitness 
v.ihies are calculated for every individual and the process is repeated. 

3.7 Non-dominated sorting genetic algorithm-II (NSGA-II) 

Multi-objective evolutionary algorithms which use non-domiixated sorting and sharing have been mainly 
criticized for their 

1. 0{mN^) computational complexity (where m is the aumbef of objectives and N is the population 



2. iinii-tilitisui approfidi. 


tli(i 1100(1 for «pe(,-ifyiiig a sharing parameter. 

•N’SCiA-II [4] alleviates all the above three difficnlties. Specifically, a fast iiou-doiiiiiiated suriiim 
ppniach with O(niN-) coinputatioiial complexity is ased. Second, a seh'clion operalur is prcseiiii'd 
liii li Croat t's a mating pool liy (rombiiiiiig the parent and child populations and st'h'ct iiig I ho bosi iwii li 
I'Siioct to fitiie.ss and .sprcjul) N solutions. In the present study NSGA-II has In'tsi used for the purpo.so 
r optimization owing to low computational requirements, elitist approach, and parametc'r-loss sharing 
pproacli. Ill the following we give a In'ief description of NSGA-II [4]. 

{.7.1 A fast uon-dominated sorting approach 

• iist . for- (’ach solution, two entities are calculated : 

1 . //,. the iiumher of solutions which dominate the solution ’i’,and 

2. .S',, a set of solutions which the .solution i dominates. 

I'lii- c.dculiU ioii of these two entities requires 0{m.N'^) comparisions. W(! identify all thosrr points which 
liavc n, --- 0 and put them in a list Fi . Fi is called the current front. Now. for each .solution in the 
|■m■l•mlt front we visit each member (j) in its set Si and reduce its nj count liy one. In doing so. if for 
;iiiy meinlK'r j the (roimt bt^comes zero, we put it in a separate list F^. When all the meinhers of the 
|■m•l•(mt. front have been cdiecked, the members in the list F\ are declared as members of the first front. 
\\'i> till'll cuiitimie this proems using the newly identified front as our current front. 

fvirli such iteration requires 0{N) computations. Since at most there can be N fronts, the worst, 
(-.■usi' complexity of this loop is 0{N^). The overall complexity of the algorithm now is 0{7nN‘^) + 0{N'^} 
or (){ in.N^). 

'.i.7.2 Density estimation 

To gel. fill estimate of the density of solutions surrounding a particular point in the population, we tixke 
ail average distance of the two points on either side of this point along each of the objectives. This 
quantity idintance serves as an estimate of the size of the largest cuboid enclosing the point i without 
including any other point in the population (we call this crowding distance). In figure 3.4, the crowding 
distance of the ith solution in its front (marked with soKd circles) is the average side-length of the the 
cuboid shown. - . ^ 



.7.3 Crowded compaxision operator 

he crowdcid (^oiiiparisiou operator (>„) guides the selection proco.ss at the various s*ages of the alyo- 
t.liiii towards a. mufornily spread out Pareto-optimal front. 
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Figure 3.4: The crowding distance calculation is as shown 

It is tissuuied that every individtial d' in the population has two attributes. 

I. .\’on-douiiuation rank (irunk) 

'1. Local crowding distance {idutancc) 

W'e now <1 (!Hu<' a partial order >„ as : 

■' J ( h'tiiik jriiiik) ***' ~ jiVrw.fc) and (idistfintT. ^ jdixtance)) 

riiat is. betwtH'u two .solutions with differing non-domination ranks we prefer the point with the 
lower rank. Otherwise, if both the points belong to the same front then we prefer the point which is 
locatc'fl in a region with le8.s(^r number of points (the size of the cuboid enclosing it i.s lai'ger). 

.3.7.4 The main loop 

Initially, a random jjarent population Po is created. The population is sorted based on non-domination. 
Fach solution is assigned a fitness equal to its non-domination level (1 is the best level). Thus, min- 
imization of fitness is assumed. Binary tournament selection, recombination, and mutation operators 
;u i' used to create a child population Qo of size N. Prom the first generation onward, the procedure is 
diff(>r<*nt. Th<? elitism procedure for f > 1 and for a particular generation is shown in the following: 

First, a combined population = Pt is formed. Then the population Rt (of size 2N) is sorted 
according to non-domination. The new p 2 urent pppulaticm Fj+i is formed by adding solutioias starting 
from the first frqnt till the sizig of the iwoiSt!^ feicofeeds the number of individdals (n) required to 



I ,)tiiiil»'t(i tli<i ])()j)uIatioii N. Tlun’oafter. the sohitious of this I'nnit an* s»trt«'<i actinli)i>; fu ami 'li« 
liisl 'll' iiidividtials aro picked. This is how we construct the population i of size .V. 'I'lim |»ojiulati<>ii 
nl' si/(i .V is now used for selection, caoss-over. mutation to create a new jxipulalion Qt . \ ul .^i/i' .V. 

3.8 Numerical integration algorithm 

A •^iiiid ordinary differential integrator should (^xert some adaptive^ control <jver its own pi’ogri’s.s. iiiakiiiv, 
li i'(iueiii. changes in its stepsiz(\ Usually the purpose of this adaptive stepsizo coutntl is to aciiieve some 
|)i<'<l<>teriuiued a,(ruracy in the solution with ininiinum computational effort. Many small steps should 
(i|)ti)e through treacherous terrain, while a few great strides shouhi speed through smooth smooth 
iniiiiter('st ing countryside. The resulting gains in efficiency aro not inen^ tens t)f pt'rcents or factors 
ol two: they can he sometimes a factor of ten, a hundred, or more. Sometimes accuracy may he de- 
manded not directly in the solution itself, but in some related conservtid quantity that can he iuonit<tr<*d. 

Im])l(nn(mta.tion of adaptive stepsize control requires that the stepping algorithm signal iulijrmation 
.ilioiit its jjerformance. most important, an estimate of its truncation error. Obviously, the calculation 
Ilf tliis information will add to the (computational overhead, but the investimuit is g(>uerally ncjiaid 
hatid.somely. 



!• igiirtc .‘{.rn Sttep doubling as a means for adaptive stepsize control in 4th (inter Runge-Kutta : points 
w here t he derivative is evaluated are shown as filled circles. The open circle represents the same 
<l('rivatives as th(i filled circle immediately above it, so the total number of evaluations per step is 11 
|)(! 1 ' two stops. 

With fourth order Runge-Kutta method, the most straightforward techni(iue by far is step doubling 
:i] as slu wu in the figure 3.8. We take each step twice, oixe as a full step, then, iud(ipendently as two 
lialf st,eps. Each of the three separate Runge-Kutta steps in the procedure requires 4 evaluations, but 
t he single and double sequences share a starting point, so the total is 11. This is to be compared not 
to 4. but to 8 (the two half steps), since — -stepsize control aside — we are achieving the accuracy of the 
smaller (half) stepsize. 

If we denc)te the exact solution for an advance from ® to x+2h by y{x+2h) and the two approximate 



snliitiDus i)o //I (c)iie Step 2 h) aud'//2 (2 steps each of size h). Since the himic nH‘lliti«l is iimrili milfi. iii« 
li ne snhil.ion and the two mimerical approximation are related l).v 

y{x + 2/1,) = y, + + 0{h.f + ... ...( I ) 

//(;,: + 2h) = {/■, + 2ihf,lj + 0{hf + ... ...(2) 

where iti ordei’ //.’ . the value (/> remains constant over the step. [Taylor series ('ximusion tells us the > ' 
is a iimuher whose order of magnitude is y('’^(;;:)/r)!.] The first equation, i.e. the (Mjuation ( 1 ) involves 
(■_'//) ’ sin<e the s1.ei)-size is 2h. while the equation (2) involves 2(/i)'’ sinc(‘ the error on each step is h ‘o 
. i’he (liff'ereiKH? between the two numerical estimates is a convenient indicator of truncation eri'or 

A = y2-yi.... ....( 3 ) 

h is this difterence that is to he kept to a desired degree of accuracy, neither two large nor two small. 
W'r ilo tins by adjusting h. 

IgiHiring terms of order h^' and higher, we can solve the two equations (1) and (2) to improve our 
iiiimerical estimate of the true solution y(x + 2/),), namely. 

y 0 r + 2 /i.) = y2 + § + O(/i)'L ...( 4 ) 

i'liis estiimit.e is accurate to fifth order, one order higher than the original Rimge-Kutta .steps. Howt'ver 
higher order is not always higher attcuracy 

.\n alternat ive step-size adjustment algorithm is based on the embedded Runge-Kutta formulas. 
origina.lly invented by Fehlberg [ 3 ]. The general form of a fifth-order Runge-Kutta formula is 

ki = hf{Xn,yn) 

hi — hf{xTi + j/m + (>21 ki ) 

.■•( 5 ) 

^6 = hf(x„ -I- aeh, pn + feet + ....+ 

Pri+l = J/n "b Cifci -h C2A2 + 03/53 -i- 04/14 -b 05/55 -f- -H 0 {h^) 

riie embedded fourth order formula is 

Vn+i =yn + clki + 4k2 + 4 kB + clki + elks + 4 kQ -b 0 {h^) ....(6) 

and so the error estmatte is . 



..{7) 


?/»+! - y;,+ , = 


rill' i)!iitic.iilar values of the various constants were found by Cash and Carp |3]. 'riiese fi;ive a iiiuie 
I'liii ii'iit method than Fehlberg's original values, with somewhat betttir error proper! ii's. 

Will'll wi' know, at least appriixiniately. what our error is, we need to keep it within desired bounds. 11' 
we i,ak(' a. step li\ and produce an error Ai, therefore, the step h,, that 'ft'ould have givi'ii .some oilier 
value is readily estimated as 

1.,, = /..^“ ...(B) 

l^ll 

1 li'iieeforth. A„ denote.s the de.sired accuracy. If Ai is lai'ger than Ao in magnitude, the etjuation tells 
us how much to decrease the step-size when we retry the present (failed) step. If Ai is smalliT than A„. 

I lieu the eipiatiou tells us how much we can safely increase the step-size for the next stej). F’urtliermori'. 
Itecaii.se our estimates of error are not exact, but only accurate to the letwling order in h. we are advisi'fl 
111 put ill a safety factor S which is a few percent smaller than unity. 

> Ai 


Ao < Ai 

This provscription is found to be a a reliable one in practice. 

3.9 GA in conjunction with numerical integration algorithm 

riie immerical integration algorithm works in conjunction with the genetic algorithm used. Initially, the 
(.I.\ which may be RGA or NSGA-II, sets a random population of design vaxiables. The limits of these 
variables is supplied to th GA code. Then every set of d^ign variables value is sent to the numerical 
integration algorithm which computes the vehicle’s dynamic response, i.e., the displacements and accel- 
erations and then transfers their maximum values back to the genetic algorithm . The genetic algorithm 
i.licn assigns a fitness after considering tha wnstrauat violations. Then the second set of design variables 





si'iit 1.0 the imuu'iiral iiitegmtion algorithm. Thus, the fitness of whole .s<>t «tf initi.'il |Mtpul;iii<tii i.s 
ili icniiiiK'd. Tlu'ii the si'loetion. cross-over and mutation opcTators op<?rate on th<? population and Idini 
a iK'w po])uhu.iou. This pui^nlation again goes through the same process. Thus, after a ninnhei of Mich 
g('neia.tious th(' optimal results are obtained. The step by step in-ocednre for th(' whole of 1 he algorithm 
is as follows : 


Step 1: Initialize the population. 

Step 2: For ('very individual of this random population, the adaptive st(ip-.size control Runge-Kntta 
method, calculates the dyuainic response for the whole period of travel for difft'rent sets of e.\(d- 
t at ion functions. The initial values of both the displacement and velocity are fixed tis zero. 

Step 3: Th(' maximum values for the displacement, acceleration and jerk are calculated. Then the 
dilfereni fr('(iu('ncies ai-e calculated, and their maximum and minimum values are found. 

Step 4: Calculate the objective function value.s as obtained from the numerical integration algorithm. 

Step 5: Rank the poi)ulation using the dominance definition defined in Sciction 3.5 rus ajiplicnble for 
iioii-dominatc'd sorting [21]. 

Step (b Calculate the crowding distance [Section 3.7.3). 

Step 7: Do selection using crowding compaiisou operator [StJction 3.7.4). 

Step 8: Do cro.ssovnr, mutation and generate the intermediate children population. 

Stop 9: C^mibine the intermediate population and parent population and do non-dominateKl sorting 
and liiid crowding distance. 

Step 10: New parent population is formed by accepting solutions from the first (best) iion-dominated 
front and continuing to other fronts successively till the new population exceeds the number of 
individuals as the parent population had. 

tep 11: The solutions from the last accepted front axe sorted according to the crowding distance and 
as many individuals are selected, which make the population of new parent population same as 
old one. 

tep 12: Repeat Step 2. 



3.10 Closure 


I'liis cli.iiiU'r iiiclii(i<\s tlu' study of various algorithms that hav(; boon tiscd iu ihc pn-soui wuik. \\t 
siail will) siini)l(' g('ucti(- algoritluu and then move on to the real-n»ded geiu'tic algoritlim. Alier siuule 
i)l)ie( tiv(' optimization using GAs we proceed to multi-objective genetic algorithms. briel'descripl iuu 
of vfuious featim^s of these algorithms is also presented. We also study an eff<H'tiv<? constraint handling 
techiiiiiiK' used in giuietic. algorithms. Then we present NSGA-II algorithm that luus bciui ns«'d in the 
present study. In the end we di.scuss the adaptive step-size control RungtsKuttii method which has been 
M>eil in eonjunetion with the N'SGA-II algorithm. 



Chapter 4 

Car suspension models 


1.1 Two dimensional model 

III :i i\v(i (Hiiii'iisioiial nuidel of a car suspension system only two wheels (one each at. rear and front i 
.ire l•ollsid('l•ed. Thus, the sprung mass is considered to have vertical and pitching motions only. .\'o 
r. illiiig iiiotion of tlii^ car can he considered in this model. Side to side .suspension interactions will he 
1 11 'gleet ed. They are of course a major factor in cornering behaviour, but the phenomena of huiincing 
iiiiil pitching are experienced in the main when travelling in a straight line. Both tlui front; and the 
re.iv wheels go though the same road ridge. The two dimensional model is shown in the figure l.l The 
iwu iliiiietisional model is better approximation than the quarter model in the .sense that tlu' pitching 
Hint inn ("in also lie considered and since the rear wheel is included, the wheellnise filtering efft'ct is also 
present in tlie responsii obtained 

1.1.1 Design variables and system constants 

III (irder t,o fonmilato the optimal design problem the first task is to identify the imijortaut design 
I ria l lies. Let ns first identify all design parameters that govern the dynamic behaviour of the car 
viliration : 

• .Si)nmg mass, m* 

• Front unsprung mass, mju 

• Rear unsprung mass, m^u 

• Front coil stiffness, kjs 

• Rear coil stiffness, krs 

• Front damper coefficient, c/ 

• Rear damper coefficient, Cr 



• Kn ml tyre slilhicss. kf, 


• Rear tyre sliffut'ss. />■,.( 

• Axle t.(i Axle distanci', 1 . 

• Momeiil of iiH'i't ia lor idlehinjj;. ./;i 


k 4 , 

Rl'iJi iinni 



l•‘i,^,•lu•e 1.1: 'File ilyiuitiile inodel of tlie car susponHion system, The atiove model has four dcffrees of 
in 'I ‘(loin 

Since a susiiensiou designer is interested in choosing the optimal dampers and suspension coils, we 
coiisidc'r only four oi‘ th<' above parameters, front coil stiffness A/.,, reai‘ coil stiffness front damper 
cocftici<!Ut cj-. rear damper coefficient Cr -- as design variables. 

.Also, most suspension coils and dampers used in cars have nonlinear characteristics. The spring 
rale or tli(> dajnping coefficient does not remain constant for all displacement or velocities. They are 
usually found to vary in a piecewise linear form as shown for the rear suspension coil in figure 4.2. For 
siinpliciiy, we consider ( when the deflection is between zero and a value (5) as the variable and other 
two spring rates, A:*, {when the deflection is more than S) and A:®, (when the deflection is negative), as 
lo vary in a fixed proportion with respect to A:®,. Dampers at front and rear are also considered to vary 
ill a similar fashion. However, the front suspension spring is assumed to vary linearly, as commonly 
f(illowt!d in automobile industries. 

Following values are obtaitaed from one automobile manufacturer : 




Krs 


Dellection 



Figure' 4.2: Biliiieai' variation of retir suspension stiffness with defliHUion 


. ].28 

• = o.oaa 

• = 0.257 


• , = 215 mill . S,.f = 50 lum/s, and 6er = 100 mm/s. 

i.1.2 Mathematical modelling 

Mndfilliiig reejiiiros several steps of mechanical and mathematical idealization. In the multi-body sys- 
ti'iii approach, deformations of the bodies are neglected and the bodies are considered rigid. They are 
(•t.mie<-t<(d by ideal links and beai’ings allowing some well defined motion of the bodies relative to (vidi 
I It her. Coupling elements like springs and dampers determine the dynamics of the systems. 


4'he ilynamic behaviour for the two dimensional model in steady state vibration can be obtained by 
\\ rit.iiig Newton’s Second law for the sprung and unsprung masses. By drawing a Free body diagram for 
each, and considering the forces acting on the sprung mass and on the front and rear unsprung mass, 
we write the differential equations governing the vertical motion of the unsprung mass at the front axle 
(i/i ). the sprung mass {q^), and the unsprung mass at the rear axle (54), and the angular motion of the 
.si)nmg mass ( 93 ) as follows (Deb, 1995) : 

91 =* {F 2 + i^3 - Fi)/mfu (4.1) 
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(j‘2 = + Fry) /ills 

<h = [{F,+Fr,)h-(Fu + F,)h]/J, 

<ii = {F\ + F, - Fv,]Jmr„ ( I- 1 i 

Tlic [)nr<un(>t<'rs l\ und /•_> nro the horizontal distances of the front and rear axli^ from llic n'lili*) nl 
,L!;i'aA-it.>- of th(' spi'niiK num. Ftn-ces F\ to F(j are caknilated as ftdlows : 

F\ = F,> — kjsd^. Fn =cfd‘^. 

F I = k‘isd\- Fry — Cy(L\. F^ — kygd'r^. 

I'lic parameters d], d-}- d;j and d.\ are the relative deformations in the front tyre, tlu' front spriiifi;. 

I lie rear tyre, and the rear spring, respectively. The parameters d‘z and ^4 are relative velocities in tln‘ 
flout and rt’ar dampers, rospectivtily. Figure 4.1 shows all the four degrees-of-freedom of the ahovi' 
svsieni (//i to f/i). The relative deformations in springs and tyres can be written as follows : 

d[ ~ q\ - flit), 

d.> =. f/._, + fqz - (/,, 

dll ~ <I\ - hit), 
d \ = q-i ~ hqii - q.i- 

The functions /i{t) and hit) are road excitations as functions of time in the front and rear tyre, 
ivspect.ively. Any suitable road profile can be tested using these functions. For example, a bump can 
111' modelled as /| (t) = A 8 in( 7 rt/T), where A is the amplitude of the bump and T the time require<l to 
cross t,he bump. When a car is moving forward, the front wheel experiences the bump first, while the 
ri'ar wheel experiences the same bump a little later, depending on the axle-to-axle distance, L, and the 
sjieed of the car V. Thus, the function hit) = hit — L/V). The sinusoidal bump profile with a width 
' if ridO mm and height 70 mm is shown in the figure 4.3 

The various system constants that are supplied as data to the algorithm are the front and rear 
iiusprung masses, the sprung mass, the itirheel to wheel distance, the distances of the center of gravity 



T 
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Figure 4.3: A sinusoidal bump profile 

tiiiiii the iVoiit and rear axle, the front and the rear tyre stiffnesses, the ear speed, the width and the 
aiiiidit.uih? of the Iniiup. These data have been obtained from a car manufacturer. 

4.1.3 Objective function 

riie two dimensional model is first optimized using a single objective real coded GA. TIk^ (jbjectiv<' 
lieii- is lo minimi/,<' the maximum transmissibility of the suspension system, i.e.. the ratio of r(vspuns(‘ 
ami)lil.ud(‘ to the excitation amplitude. The response anrplitude in this case is for both pitching and 
honnc(' motion. This optimization is done subject to only one constraint, i.e., jerk constraint. A v(‘ry 
im|)<»rta.nf. featun? of this optimization is that we are formulating the objective function in two ways, in 
t lie first cas(f we minimize the transmissibility only for the period when the car is not on the hump, i.e.. 
cxclnde the time when any one wheel of the car is on the bump. In the second ca,se the transmissibility 
i.'i minimized for whole of the period. 

However, as single objective optimization considers only one objective from a set of objectives, and 
o\erlooks the others, the optimal results obtained are very often not too good from a practical point 
uf vi(!w. The next step of optimization is multi-objective optimization. Two conflicting objectives, i.e., 
minimizing tin? maximum transmissibility and minimizing the maximum acceleration are considered 
rsnh ji'ct to only one constraint, i.e., jerk constraint. The optinaization run provides us with a set of 
i optimal solutions with a reasonable trade oflf. The designer can choose any of these solutions according 
'to specific requirements. 



1.1.4 Constraints and variable bounds 


III till' r.wi) (liiiu'usii)ual model only one constraint is consuleriHi. i.o.. jerk constraint. 'I’lic optimal <li’- 
si 5 >;ii must not violate this constraint. The constniint speiafies the jerk (the rate of chalice of veriirai 
accelerat ion) t.o lie le.ss than 18 m/s'* (as suggested by an automobile manufticiurer). The viTlical Jerk 
can 1 ) 1 ' comiiuted by numerically differentiating the vertical acceleration of th<‘ sprtiug mass 

Aiioi lier important constraint is the restriction on the rattle space. It would be nice to have as soft a 
sitspiMisiou as ])ossible. That means lower stiffness values but this may lead to crushing of the car parts 
liei wemi the axle and the chasis. So there must be a minimum rattle space constraint. .Also, a softer 
suspension leads to problem in the handling of car, so from the safety point of view it is imperative that 
I lie snsiieusion is not too soft. However, these things can also be handled by putting a lower botind on 
I lie stiHiiess and ilamping coefficients. So, we have used some lower and upper bounds for the <lesign 
variables. 

-1.1.5 Algorithm used 

['’or t he multi-objective optimization, a multi-objective algorithm, by the name of Non dominated sort ing 
genet ic algorithm -II (NSGA-II) [4] is used in conjunction with a numerical integration technique. 'I'lie j 
coupled diffen'iitial e<iuations specified in equations (4.1) to (4.4) are solved using adaptive step-size I 
control Rimge-Kutta fourth order algorithm. The numerical integration routine finds the dynamic i 
response of tin* susi>ension system to the road excitation for a set of design variable values provided by | 
I he genetic algorithm and returns back the maximum transrnissibility and maximum acceleration values | 
which become the two objective's for the multi-objective algorithm. ‘ 

In the (tase of single objective optimization, a real coded GA with Simulated binary cross-over is 
used ill conjunction with the numerical integration technique. The real coded GA works in t.iie real 
space so more accurate results are obtained. 

4.2 Three dimensional model 

In t he t hree dimensional model, all four wheels are considered. Thus, the rolling motion of the sprung 
mass can also be studied. The sprung mass can have three motions — vertical bouncing, pitching, and 
lolling. Btisides, each of the four unsprung masses will have a vertical motion. Thus, there are a total 
of seven second-order differential equations governing the motion of the sprung and the unsprung masses. 



•1.2.1 Design variables and system constants 

Siiict' both U>f't and right rear (or front) wheels have the same suspension system, flie iminber ol doiuii 
variables in the thret! dimensional optimal design model is also four. However, the «lynanii<‘s ul iIk 
ear will lx- dilhnent from that in the two dimensional case. Also, there are two mon' system eonslaiils 
besides I h(' systenn eonstants that are present in the two dinumsional model. uauK’ly: 

• MdiiK'iit of inertia for rolling. 

• Dist,anc(^ Ix^tweon th(> loft and right wheels, Z,. 

In ;i-D 1110(1(^1 also, the l)ilinearity of the dampers and the niar springs iwe considered. They vary in 
,1 pieeewi.se liiu'ar fashion similar to that in the 2-D model. However, we consider the luodiKed valmw ' 
of the left and the righti suspension separately as the displa(;ement and velocity values for the two vary i 
dilCerently. Tlx’ limits on the value of deformation and the rate of deformation in the springs and t he 
(lain])ers respectively are the same as used in the 2-D model. 

1.2.2 Mathematical modelling 



Figure 4.4: A dynamic car suspension model with seven degrees of freedom 







ill tilt' ;J-D uiDtli’l as sliown in figure 4.4 there arc a total of seven st*eon<l onler ililfereinifil equal it ni> 
whitii are derivetl using Newton's second law for the sprung and the unsprung in;is.ses. M.v tlrawiiiu a 
frtiti I>t)tly diagrani fur each, and considering the forces acting on the .sprung niiuss anil on the li'uiit ami 
ifar unsprung mass, we write the differential equations, for the vertical motion of the uns|)rung masse*, 
,il tlie rigiit siile ami the left side on the front fixle (</| and t/o respectively), the vertical motion of tin- 
sprung mass the vertical motion of the unsprung in<i.s.ses at the right side ;ind the left siiie on the 
mar .Lsie (tjr, aiifl qj iesj)cctively) and the pitching and rolling motions of tiie .sprung nutss (t/.j ami //| 


l•l‘Sl)l't•t,iv('ly). The' equations are as follows ; 

(ji - (F-i + F^- Fi)lm.iu ( I-'*) 

<h — + Fz + F.i + Fr, + Fa + Fi) + Fiq + F\ \ (4.0) 

<h = [(■P’t + -P 5 + + Fi\)l2 — {F 2 + Fi + Fg + Fa)li]/Ji, (-1.7) 

= [(P's + -Pg + ■P’lo + Fii) — (F 2 + + Fi + F 5 )]/, (4..S) 

95 = (F} + F5 - Fa)/mru (4.9) 

9 i) = (Fg + Fg - F-!)lin[u (4.10) 

f /7 = (Fio + Fii - Fi2)/rnru (4.1 1 ) 


rite forces F| to F12 are calculated as follows: 

F| = kftd\r Fa = kf^da, F3 = c/tia, F4 = A:,.sd 4 . 

P'r, = c, .<•/.(. Ffi = krtdz, Fr — kjtdz^ Fg = kf^d^, 

Fg = C/dg, Fio = krads, Pn = Crdsi P12 = krtdy. 

The paxaineters di, da, dz and d4 are the relative deformations in the front right tyre, front right 
spring, rear right tyre and the rear ri^t spring respectively. The parameters dg, de, dr, dg ar® rela- 
tive deformations in the front left tyre, front left spring, rear left tyre and rear left spring respectively. 



1 'll (/,.>. </). f/i, aiul tls art* rclatiw vtilocities iu the front, right, (laiiipcr. rear right (laiiiper. 

Ii, .iii ii‘i! ilaiiipfi ami rear Iclt <laiupi'r respect iv<*ly. Th<^ relativ<^ defonuations in th*' springs ami the 
i\ie> i-.ni lie written as folhiws ; 




./ - /i</a *■ /;(/)- </i. 


'/.! - -- hV). 

(/ 1 f/j ■ Ijt/.i + /,.«/] - (fr,. 
'/•. - '/I. ,/l(/). 


du ~ II , + /,//;$ - /,.y.i - //,j. 


d: '/: n. 

f/s - q, /jf/a - 

Tlie fiinctions f\(t) ;uu! /^{t} are road excitations as functions of time in the front and rear tyre. 
i espect iv<*ly. ’Pile roaii proHU? that was tested in the two dimensional case was a sinusoidal profile. In 
' t hree dimensionnl ctust* different excitation functions used are as follows : 


I. A hmnp motlelled as fi(t) = di[l-co8(27rt/r)]/2, where A is the amplitude of the bump and 
T t.h<' time nxpiired to cross the bump. The rear wheel experiences the same bump a littU' 
biter, dtipeuding on the wheel to wheel distance and the car speed. Thus, the function f^it) = 
!\{i - LjV). A bump of this profile having a width of 2000 mm and height of 70 mm is shown in 
figure 4.5. It is to be noted that the same kind of excitation is experienced by the right and the; 
left tyres and we ignore the lateral interaction between the right and left suspensions. 


1. Three similar bumps dually sptffcced havuig tfie sme profile as in above case. The gap between 
the bumps is taken as ^ual to the tridth erf the btimp. The bumps in series are shown in the 
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Figure 4.5: Single bump profile 


ligiire i.H. 



Figure 4.6: Three bumps in series 
A rounder! displacement step function as follows [15]: 

/){/,) = A[l — (1 + 'yT)e^~'^^] where A is the amplitude of the displacement step and T is the 
uon-dimensionalised time, i.e., T = t/{L/V) and 7 is the severity parameter for shock excitation. 
This kind of base excitation with 7 equal to 0.1 and 50 is shown in the figure 4.7. This type of step 
function is used to see the effect of noa-linear damping. For the rolling motion, we have considered 
a rounded pot-hole which has the same mathematirial expression as above, but with a negative sign. 
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Ki>;nr(! 4.7: RouiukKi displacement step function 

1.2.3 ObjtH-tivc? function 

hirst ly. i'cm! cudiMl (JA is tiscd to optimize a single objective, which i.s considentd to be the vertical 
ir.iiisiiiissibiiiiv. in the <-ase of vertical botmce motion. 

Ill the fa.se of pit«-hitig and rolling motions the single objective considered is pitciiing and rolling 
.(iniilii tides of the sprung innss respectively. In the case of 3-D model we consider only the case in 
uliii'lt the whole lime of travel is cousidttred. as in the case of transient response it is moiv importiint, 
to niininii/.e ilie maxitnmn values of displacement ami acceleration even if it comes at the cost of little 
ride romlnrl . 

Aitei the single objective optimization multi-objective optimization is performed by considering two 
foiilih ting obj<*etives which are of importance to the suspension designer. The sutspension .system is 
siil.jeeted to all the above four set of excitation functions to study the performance of the .suspension 
l■llilrilcleristi(■s to the different road profiles and to understand, the control that has to be exoii.ised 
nil t lie difiereiit desigti variables in order to obtain a suspension system that is good for a wide range 
nf road <>xcitHtions. The various conflicting objectives that have been considered for the purpose of 
imilti-obj<!ctive optimization are as follows : 

1. Th(‘ vertical bounce motion and the vertical acceleration. We have considered the vertical trans- 
missihility in place of vertical displacement. The transmissibility is given by ^2 

2. The pitching amplitude and the pitchii^ acceleration. 

3. The rolling ampUtude and the rolling teC^eratibn. Hbttrever, the tollihg motion has been Consid- 
ered in the case of the rounded displacement step kind of fiinction only. B is used to represent 



Mill' M'l 111 wIh'i'Is Iriuln 


or li'll ) tailing in u pot-hol(» which arc so comnioiily tonud on Ihc mads. 


1.2.4 C’oiistraints 


111 ilic iwti diinciisioii.d model oulv tin* jerk constraint wjvs considered, 
tiinv nv live lonsiiaints which are considered, 'i’he optimal design 
I Miisiiainis. 1 lie et msi I ainl.s are eniinierated Ixdow : 


In the three dimensional imxlcl 
nmst not violate any of these 


I he verlie.'il jerk inns! not exceed IK m/s'*. The vertical jerk is calculated by numerically difl'er- 
I'liiiaiiti'j. the vertical ac«'eler;ttion of tlu‘ sprung mass c/a- Mathematically, it is uxpres.s<Hi as: 




tUfj 

li 


< 18 


This const t aint is inelnd(‘d for tlu* rhle comfort consideration. 


•J. The nia.vimnni of the front, rear and pitching frequencies should not exceed 2.5 Ez. The higher 
the lit'<piencie.s the iintre uncomfortable the ride. The lower accelerations occur at lower natural 
lietineiieies .So the ciir suspension should not have very high vahuis of frequencies fis tluty Ictad to 
highei ai i-elei atitm valims. Mathtunatically, it can be expressed as 


Fw < 2.5 

:i. riie ininiimim of t lie front, rear, and the pitcliing frtxpiencies shotild be greater than 0.8 Hz. The 
lower lieipieneies lead to mort? ride comfort, but they are obtained by using softer suspensions 
which are not .so good from the ride control point of view. So a lower bound is fixed for the 
frecptem'ies to htive better ride control.Mathematically, it is expressed as : 


1. .Also I he front natural friniuency should be at least 10 percent less than the rear natural frequency. 
This is dotie to itiakt* the pitching motion die faster. This may lead to longer bounce motion but 
vertical bounce is lt«is annoying as compared to the pitching motion. 

fe®, > l.lKf, 

5. Another important constraint is the constraint on the maximum value of acceleration. This is 
(^specially important in the case of muH^bJectiVe optimization as in the optimal results certain 
acceleration values exceed 9.8 m/s* (which is the immieiricai -value of aecelCTation due to gravity, 
i.e. 'sf). Since Wt sudi acceleration values the passenger will lose contact of his ^at, snch values 



.uf ii.M'Ii >>■ III ill'' |>ri'M'iii sluiiy a const raiu< restrict injt; the niaxiinmn .icci'lcratioii value to !).li 
111 >■ i' HM « 1 . Mat liemat ically. the const raint is written as : 

< !).0 

lli^iilc", ihi'N' (unsliaints. the variai>ie hoinuis a,re also fixed in order to atx'ount fur oilier coii- 
^naini'. 1 he niiniiiiiini rattle space* constraint is considered by putting a lower hound on the 
■.iiirnes> values. This also lakes can* of the better ride handling aspects which are a must for the 
j ii ir V 

1.2..’) Noii-linojirity in damping parameters 

1,; the 0-1 ) model. I he dampers with iMibic uon-liuear damping tu’e considered [C]. when the road prolile 
i- iiiodelled .<> a lonuded dispiaceiiieut step. The linear model is ustially considered for mathematical 
- iiiiiliciiy Imi ill «‘ase iln* transient tlisplaceunnits are sufficiently large, then the non-linearity in the 
^;llill•J^ utd the daiiipers caimot ht* ignon'd. The non-linearity in the dfunpiiig rather than that in 
-'iilii. ■ss has a mole pi ill toll no *ii effect on dill'ereut iudietjs expressing the performance of a suspension 
-vsieiii. So. ill the present study we imve cottsidered cubic non-linear dmnping in addition to linear 
ilaiiipiiig. while only linear stiffness characteristics are considered throughout. In the present work we 
h.ive taken a negative ilauiping coeHicieut which has a numerical value equal to oiui-fonrth of the value 
III lineal damping coellicient. 

1.2.() Matlicmatical formulae used 

The pitehing and t he rolling moments of inertia are calculated as follows: 

• I 'olar nioinenl of inertia 

Jp = m/ii + mrll 

• Rolling moment of inertia 

j; = (m/ + mr)Ir 

For t ill? abtive formulae of mooients of jinertia, m/ and ttv i*!''®! front and rear sprung masses. 



[ ill 1 .ii< iii.ii i<*u "f lilt' Iroiii anti Jln> it*av natural tmjuoudos tiupeiuis on I.Ik' center ol' gravity ami 
jiiiii'j. leiiit'i "i lilt till ami Ollier parameters of the ear suspension system and is derived elsewliere 
. I Ilf Iroiii I / ■ ! and rear i h, ) natural frecjueney and the pitching natural frecjueney (F^,) are given 
.. f, .ilo\\> : 


i l‘ 2 kj,{a + .s - x)-^ + 2 A:<;,(.s - /) - .r)- 
2;ry“ ' ■ /n,,(/v2 + .s^') 


I l'2k/siu - r - x)‘^ + 2Ar{'^,(r + b + x)^ 
2 “ y + 7 -^) 


r L ~ + x)'^ 

''' 27r V m,K‘^ 


where I is the ilistauce hetwisut the center of gravity and the spring center, o and h are the dis- 
i.iiices of the spring center IVtuii front and rear wheels : 


a = 


Jfc® 

—HUL—Ij 
^rs + kfs 


6 = 


,—hh—lj 

kfH + kfs 


The paranieters r and s are the distances of the front and the rear conjugate points from the spring 


n'uH'i ^iv(‘n as Inlhiws : 


K^-x?-oh 


r = 


2 x 



and K is tlin radius of gyration about the center of gravity. 

4.2.7 Nondimensionalization of the differential equations 

Siu(;(! w(' are using the adaptive step-size control Runge-Kutta method, we must guard against the stop- 
si/<' Ixsauuing too small, otherwise the routine will make the apparent error zero and chug happily along 
taking iufiuitoly many steps and never changing the dependent variables one iota. This may sijecially 
liai)i)en at higher values of speed. Non-dimensionalisation of the equation helps to avoid this. In the 
i lire(^ dimensional case, the differential equations were non-dimensionalised. For illustration a few (jf 
I lu! ('(luations are written below and their non-dimensionalised forms are also written alongwith. All 
till- non-dimensional terms are written as capital letters. 

T = — ^ 

L/V' 

^ _ h{t) 

A A A ' 

Di = Qi-Fiit), 


_ g2 . h<l3 . IrQA ^ 
A~ A^ A A “ a ’ 

r>2 = <32 + + ^<34 - <3l» 


^2 _ 92 , hqs , ir94 qi 

A A'^ A ^ A A' 

D2 = Qi + ^-jQz + —<34 - <3i, 
A a 


An example of the non-dimensionalised form of differential equation is given alongwith the original 
differential equation. Again the non-dimensionalised terms are in capital letters. 

q'l = {kfsda + c/da - kftd\)lmfy,, 
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Qi 


fc/,L2 

F2m/„ 


D2 + 


CfL 

Vm/u 


D2- 


mfuV^ 


Di 


'I’liis way all the clifierential equations are non- dimensionalised ajid used in the numerical integra- 
I ii)ii routine. The non-dimensionalisation helps in proper working of the routine even for higher valiuis 
I )!' the cfu' speed. 

4.3 Closure 

'I'liis cliaptoi- includes the discussion on two and three dimensional models. Mathematical modelling is 
presented and the problem formulation is presented for both the cases, i.e. the objective function and 
the constraints are defined. In the three dimensional case, a discussion on nonlinearity of dampers has 
l)(sni included. The non-dimensionalisation of equations in the 3-D case has been discussed in the end. 
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Chapter 5 

Results and Discussion 


Till' siiigle-()hj(u;tivfi optimization and the multi-objective optimization ai'e performed on the two and 
(linu'nsional models and the results obtained are qualitatively studied. 

5.1 Two dimensional model 

’riie following pajameter values are used in all simulations; 

III, - 7.‘3() kg. rrifu = 50 kg, 


111,11. = 115 kg. 


kft = 15 kg/mm, 


A',, = 17 kg/mm, 


h = 1.50 m, 


/.. = l.iin m. 


L = 2.85 m, 


I ' = 5 kin/hr, 


Jp = 2.89(10®) kg.mm^ 


The (;ar motion is simulated over a sinusoidal bump having 500 mm width and 70 mm height. In 
all .solutions, the spring rates are expressed in kg/mm and damping coefficients are in kg.s/nun. Tlie 
vai iahk' bounds for the stiffness values are 0.50 to 8.0, while those for the damping coefficients are 0.10 
to 8.0. We start by presenting the single objective (RGA) results for the two dimensional model and 
t hen follow it by the results of multi-objective optimization. 

5.1.1 Single objective optimization results 

The real-(»ded GA is used for single objective optimization. The various parameter values are as follows ; 
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Pi>piilation size : 100 
(.'i-oss-over probability : 0.85 
Mutation probability : 0.25 
.\-over (listriluitioii index : 20 
Mutation distribution index : 100 
•N’ umber of generations ; 200 


Results with bump time excluded 

111 the singl(! objective optimization results we have considered two cases in which the objective function 
is evaluated differently. In first case the objective function i.e. the vertical transmissibility is considered 
only for that part of cai' travel when none of the wheels is on the bump. In second case, we consicler 
the wliole time of travel, for the evaluation of the objective function. In both the cases the objective 
limctioii is subjected to only one constraint, i.e. the jerk constraint. We will first present the results 
liu' tlu' former case. 




figure 5.1: Coinparision of present and previous results for the case in which the bump time is excluded. 


Results 

kfs 



4 

Transmissibility 

Present 

5.73 

0.73 

3.26 

0.50 

0.66 

Previous 

(TELCO) 

1.56 



1.00 

0.77 


I'able 5.1: Comparision of optimal values of stiffness and damping parameters for the case in which 
bump time is included. 
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Thf^ results niai’ked as present are obtained in the present study , while the previous results are that 
used by an automobile manufacturer. It can be clearly seen from figure 5.1 that the present results are 
l)ett.<'r in terms of both the displacement and acceleration. The optimum values of stiffness mid damiiing 
l)ara.nieters me given in the table 5.1. The maximum vertical jerk as calculated for the present study is 
17.9 ni/s'^ which makes the jerk constraint an active constraint, and indicates that the obtained result 
is the global optimum. 

Also, we will compai-e the results obtained in the present work (using RGA) with those obtained 
l)y using binary coded GA (SGA) which are presented in M.Tech thesis of Mr. Vikas Saxena [1]. A 
ciunjiarision of results with binary coded GA is presented below. 




Figure 5.2: Comparision of real-coded GA and binary-coded GA results 


Results 

kfs 


ua 

^TS 

c? 

Transmissibility 

RGA 

5.73 

0.73 

3.26 

0.50 

0.66 

SGA [1] 

4.53 

1.72 

2.86 

1.01 

0.71 


Table 5.2; Comparision of optimal values of stiffness and damping parameters obtained by using RGA 
and SGA . 

It can be seen from the figure 5.2 that the results obtained using RGA axe slightly better than those 
obtained previously by using simple genetic algorithm with binary coding. Both the displacement and 
acceleration responses are better in the present study. The values of stiffness and damping parameters 
for both the cases is given in the table 5.2. 
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Results with bump time included 

Results of second case, in which the objective is m i nim i z ation of the maximum vertical displacement 
lor tli(i whole time of travel, i.e, even the time when the wheels are on the bump ai'e considered, ai-e 
now presented. We again compare the results with those used by one of the automobile mauufacturcu' 
(TKLCO). Even in this case the results obtained are better as can be seen from figure 5.3. In table 
th(^ values of stiffness and damping parameters for the two cases is presented. 




Figure. 5.3: Comparision of present and previous results for the case in which whole of time of travel is 
considered. 


Results 





Transmissibility 

Present 





0.77 

Previous 

(TELCO) 

1.56 



1,00 

0.93 


'Table 5.3: Comparision of optimal values of stiffness and damping parameters for the case in which 
whole of time of travel is considered. 

An interesting observation can be made when the results of the first and the second case ai-e com- 
parcMl. It is observed that if the vehicle suspension is designed firom the point of view of failure under 
stress, the second case is better but the first case leads to faster damping of vibrations. A comparision 
of the two results is presented now in figure 5.4 

Looking at the values in the table 5.4 it can be seen that for the case in which bump time is 
('xcluded we have larger stiffness values as compared to those in the case when the bump time is also 
included while considering the objective function. So, in the former case we have higher displacement 
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5.4: Comparisiori of two cases: one with bump time excluded and the other with bump time 
iiiclncic'd. 


Ill- t.raiismissibility values. Since, the damping is higher in the first case we have faster damping in the 
first case. 


Results 

kfs 




Traxismissibility 

Bump time excluded 

5.73 

0.73 

3.26 

0.50 

0.91 

Bump time included 

1.23 

0.36 

3.28 

0.24 

0.77 


I’ahle 5.4: Comparision of optimal -values of stiffness and damping pareimeters for the two cases in which 
Immp time is excluded or included. 

5.1.2 Multi-objective optimization results 

III th(' present study we have considered the minimization of transmissibility in bounce motion and 
iiiiiiimization of vertical acceleration as the two objectives. The results were obtained in a single run of 
t Ik' NSGA-II algorithm used in conjunction with the adaptive step-size control Runge-Kutta algorithm. 
The various parameter values axe as follows : 


Population size : 300 
C'ross-over probability : 0.85 
Mutation probability : 0.25 
X-over distribution index : 20 
Mutation distribution index : 100 
Number of generations : 200 
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Results with car speed 5 km/h 


111 the multi-objective optimization we consider the case in which the bump time is also included, i.e. 
f.lui whole time of travel is considered while evaluating the objective function. The reason for this is 
that in the transient vehicle response we are trying to minimize the maximum values of acceleration 
and displacement as they are the critical factors in the failure of parts under stress. The results of 
multi-objective optimization are presented now. We first show the set of optimal values in the objective 
function space in figure 5.5. These values are called the Pareto-front. 



Transmissibility 


Figure 5.5: The Pareto front showing the optimal valu«i in objective space for the car speed of 5km/h 


Now, the designer can choose any set of values from this set depending on his specific requirements. 
The range of values can be better observed by looking at the two extreme values. The response based 
on these two extreme values is shown in the figure 5.6. 


Results 

kfs 

m 

H 

■ 

Transmissibility 

Max acceleration 
(m/s2) 

Minimum transmissibility 

0.50 

0.10 

0.69 

0.327 

0.60 

2.8 

Maximum transmissibility 

0.53 

0.13 

0.669 

0.135 

0.84 

1.27 


Table 5.5: Optimal values of stiffness and damping parameters with car speed 5 km/h 


Prom table 5.5 we can see that the range of acceleration and transmissibility is not large in case of 
car speed of 5 km/h. The reason is that we have used a very low speed, i.e. 5 km/h and at such low 
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5.6: CoiiiparLsion of two extreme values 
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Time of travel (s) 

the Pareto-front for the cai’ spee<l of 5 kiii/h 


spiKids there are not many variations possible and mostly it is one solution instead of a set of solutions. 
Results for single bump excitation at a car speed of 25 km/h 

I'Voni figure 5.7. it can be seen that if we change the car speed to 25 km/h, other parameters remaining 
tli(' same a wider range of solutions can be obtained as compared to those at car speed of 5 km/h. 
I’igun* 5.8, clearly shows the comparision of two extreme values of the Pareto-optimal front of the 
iigure 5.7. In table 5.6, the respective values of stiffness and damping are given. 


Results 

kfs 



c? 

Transmissibility 

Max acceleration 
(m/s2) 

Minimum transmissibility 

1.86 

0.11 

2.05 

0.23 

0.70 

3.52 

M aximum transmissibility 

0.62 

0.14 

0.52 

0.10 

1.12 

1.40 


Table 5.6: Optimal values of stiffness and damping parameters with car speed 25 km/h 


5.2 Results for three dimensional model 

1)1 the three dimensional model the objective function is subjected to three more constraints, i.e fre- 
(luency constraints as was discussed in chapter 4. Also, when the vehicle speed is 50 km/h, we add 
iiuother constraint for limiting maximum acceleration. The left and right suspension parameters are 
considered to be the same. The car parameters used in this model are the same as that used in the 
t wo dimensional model except = 1460 kg, the left to right wheel distance is 1.462 m, and the rolling 
moment of inertia is Jr = 1.864(10® )kg.mm^. Also, the width of the bump is taken as 2000 mm. 
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Figure 5.7; The Pareto front showing the optimal valu^ in objective space with a car speed of 25 kin/h 
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Figure 5.8: Compaxision of two extreme values on the Pareto-front for the car speed of 25 km/h 


5.2.1 Single objective optimization 

In three dimensional model, different kinds of road excitations are considered. In the case of three 
dimensional model, we have also considered the pitching and the rolling motion when the car is excited 
by a rounded displacement step. Also, the profile of the bump has been dianged from sinusoidal to 
profile represented by 1/2(1 -cos(27ra;)) Also, for all the cases two cax speeds are considered, 5 km/h 
and 50 km/h. 
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Results for single bump excitation at a car speed of 5 km/h and 50 km/h 

iMrut.ly. the results for the case when the car experiences a single bump of the profile 1/2(1 -cos(27r;r)) 
arc i)rcs(Mite(l in figure 5.9. 



Time of travel (s) Time of travel ( s) 

Figure 5.9: Results for single bump excitation at a speed 5 km/h 

When the car passes over a bump of width 2 m, with a speed of 5 km/h it is observed that the 
s])nmg mass almost follows the bump, and the wheelbase filtering effect is not very predominant, as by 
the time the rear wheels come in contact with the bump, the vibrations due to the first wheel hitting 
I he bump have almost died down. But, when the car speed is 60 km/h, the wheelbase filtering effect 
t:fui be clearly seen. The vehicle response for single bump excitation in the case of car speed of 50 km/h. 
is shown in the figure 5.10. 


Results 

kfs 

cj 

La 

^rs 

c? 

Transmissibility 

Max acceleration 
(m/s2) 

Speed 5 km/h 

1.17 

2.76 

1.24 

5.21 

0.56 

0.62 

Speed 50 km/h 

1.49 

1.05 

0.86 

0.26 

1.48 

9.00 


Table 5.7: Optimal values of stiffness and damping parameters for single bump excitation 

By studying table 5.7, we see that in the case of car speed the damping coefiicient values are much 
smaller in comparision to those at car speed of 5 km/h. This can be understood because at higher 
speeds the higher the damping force the greater the disturbance fed into the sprung mass by any given 
load irregularity. For the speed of 5 km/h with the single bump the values of front, rear and pitch 
natural frequencies are 1.21 Hz, 1.35 Hz and 1.29 Hz respectively. Also, the maximum jerk value in 
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Figure 5.10: Results for single bump excitation at a speed 50 km/h 

t his case is about 6.4m/s^. In the case of single bump excitation at car speed of 50 km/h the values 
1)1' maximum vertical jerk, front frequency, rear frequency, and pitch frequency axe 14.9 m/s’^ . 1.13 Hz. 
1.37 Hz and 1.25 Hz respectively. So, in both the cases there is no constraint violation. 

Results for three bumps excitation at car speeds of 5 km/h and 50 km/h 

Next we consider three similar bumps in series, with the gap between the two bumps same as the bump 
width, i.e. 2000 mm. Again the optimization is performed for two speeds 5 km/h and 50 km/h. 







Looking at the figures 5.11 and 5.12 for the three bump excitations at two speeds, we observe 
t liat for speed of 5 km/h there is little wheelbase filtering while in the case of speed 50 km/h there is 
sigiiificHiit wheelbase filtering. Also, at higher speed the displacement and acceleration values are highei’ 
.■IS ('xpected. 




Figure 5.12: Results for three bumps excitation at a speed 50 km/h 

Now looking at the values of stiffness and damping parameters in table 5.8, we see that at higher 
sp<‘(‘(l the stiffness values have increased and the damping -values have reduced considerably. It is b(> 
cause the higher stiffness values are required to reduce the displacement and lower damping values do 
not feed mudi disturbance to the sprung mass. The values of front, rear and pitch natural frequencies 
in the case of car speed of 5 km/h are 0.84 Hz, 1.35 Hz and 1.14 Hz respectively. The values of front, 
rear and pitch natural frequenci^ for car speed of 50 km/h are 2.08 Hz, 2.38 Hz, 2.22 Hz respectively. 
Since, lower frequencies correspond to higher ride comfort it is obvious that the ride is more comfortable, 
when the car speed is 5 km/h. However, to achieve such low frequencies, a softer suspension has been 
lusiid, which asks for a compromise in ride handling, but then ride handling is not a problem at low speeds. 


Results for rounded displacement step excitation at car speed of 5 km/h and 50 km/h with 
linear and non-linear damping 

To study the effect of non-lineax damping a more severe excitation function is required. So we have 
considered a rounded step function and a cubic damping term is added to the linear damping term in 
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■ 

Transmissibility 

Max acceleration 
(m/.s2) 





0.86 

0.54 

1.60 





0.68 

1.10 

6.50 


Table 5.8; Optimal values of stiffness and damping parameters for three bumps excitation 

I li(! force; expression. The coefficient of cubic damping term is taken to be negative and one-fourth of 
tlu; lineal’ damping term. We now present a comparision of the vehicle response with linear damping 
and cubic damping at two speeds 5 km/h and 50 km/h. A negative damping coefficient with numerictal 
valiK* one-fourth of that of linear damping coefficient has been employed. 




Figure 5.13: Results for rounded step excitation at a speed 5 km/h 


It is evident from the figures 5.13 and 5.14 that the non-linear damping with the negative co- 
(;fficient has worsened the -vehicle’s response. This non-linearity effect is more dominant at higher speeds. 
The values of stiffness and damping parameters for linear and cubic damping for the car speeds of 5 
km/h and 50fcm/h are presented in the tables 5.9 and 5.10 respectively. 


5,2.2 Results of multi-objective optimization 

.lust like the two-dimensional case we consider the two conflicting objectives the transmissibility in the 
liounce motion and the vertical acceleration. Also, other conflicting objectives like pitching amplitude 
iiud pitching acceleration and rolling amplitude and rolling acceleration are considered. In all the cases 
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Results 



jua 

^rs 

c? 

Transmissibility 

Linear damping 

2.35 

0.47 

3.72 

0.98 

1.34 

Cubic damping 

1.41 

0.59 

3.82 

0.95 

1.40 


l;i.l)I<; 5.9: Optimal values of stiffness and damping parameters for rounded step excitation at speed 5 
kiii/li 




Figure 5.14: Results for rounded step excitation at a speed 50 kin/h 
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kfs 


Ks 

< 

Transmissibility 

Linear damping 

4.77 

0.25 

3.02 

1.50 

1.74 

Cubic damping 

3.68 

0.29 

1.14 

1.59 

1.92 


Table 5.10: Optimal values of stiffness and damping parameters for rounded step excitation at speed 
50 kni/h 


the whole of the time of travel is considered. We observed that when the car speed is 5 km/h the range 
of transmissibility values is very small, in the case of single bump and a series of bump. This happens 
biM-ause the car speed is very small and we have increased the bump width to 2000 mm, so the car 
follows the bump smoothly and not many variations are available. So for this speed the RGA results 
n,r(* the best. 

Results for single bump excitation at a car speed of 50 km/h 

Now, we present the results of multi-objective optimization when the car is excited by a single bump 
which is of the same profile as was used in the single objective optimization. The speed of the car is 
50 km/h. We start by presenting the Pareto-front obtsdned in this case in the figure 5.15. The two 
{^xtreme solutions are plotted in the figure 5.17. 
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Figure 5.15: The Pareto front showing the optimal values in objective space for single bump at 50 kin/li 

Figure 5.16 clearly shows the exact position of the Pareto-optimal front in the feasible search space. 
The small dots in the figure represent the points which do not violate any of the constraints. The 
imdti-objective optimization results presented in the table 5.11 point to a very interesting fact about 
the damping coefficients. It is seen that when acceleration is considered as one of the objectives the 
damping coefficients tend to move to the lowest minimum possible most of the times. The reason prob- 
ably is that higher damping coefficients that the higher the damping force the greater the disturbance 
fed into the sprung mass by any given road irregularity, so the optimal solutions have damping values 
which are just enough to prevent undue persistence of vibrations. 


The results presented in the table 5.11 clearly show that the higher stiffness values lead to lower dis- 
placements and higher accelerations and a softer suspension produces smaller accelerations. The value 
of front, rear and pitch frequencies corresponding to noinimum transmissibility are 1.44 Hz, 2.50 Hz and 
1.99 Hz respectively. So, in this case the maximum frequency constrmnt is active assturing us that the 
results obtained have a high probability of being the global optima. The front, rear and pitch frequency 
values for the maximum transmissibility case are 0.80 Hz, 1.57 Hz, 1.22 Hz respectively. Again, in this 
case the minimum frequency constraint is active and thus it can be said that the results are optimum. 
Also, the stiffer suspension corresponds to higher frequency values and the softer suspension corresponds 
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Figure 5.16: The position of the Pareto-optimal front in the feasible search space in case of single bump 
excitation at speed 50 km/h 




Figure 5.17: Results for single bump excitation at a speed 50 km/h 
to lower frequency values. 
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Results 

kfs 


jua 

^rs 

■ 

Transmissibility 

Max acceleration 
(m/s2) 

Minimum transmissibility 

4.95 

0.13 

1.40 

0.10 

0.58 

6.3 

Maximum transmissibility 

1.97 

0.17 

0.44 

0.16 

1.19 

3.75 


'rjU)l(! 5.11: Optimal values of stiffness and damping parameters for single bump excitation at speed 50 
km/h 


Results for three bumps excitation at a car speed of 50 km/h 

Next, we consider the excitation of car suspension with a series of three bumps equally spaced with the 
gap between them equal to the width of the bump. The car speed is taken to be 50 km/h. 

The Pai-eto-front is shown in the figure 5.18. Also, the position of the Pareto-optimal front in the 
search space is shown in the figure 5.19. 



Figure 5.18: The Pareto front showing the optimal values in objective space for three bumps in series 
at 50 km/h 

The two extreme responses obtained from the Pareto-optimal front are shown in the figiue 5.20. In 
case of excitation of car suspension by a series of bump, the vibrations take a long time to die down, in 
the case of lower transmissibility and higher acceleration. The reason is that the vibrations produced 
by a single bump are further accentuated by other bumps in series. But, for higher transmissibility 
case, the effect of other bumps does not add to that of the first bump as the sprung mass has already 
been disturbed to a very displacement value. The two extreme solutions are given in the table 5.12. 
For the case of minimum transmissibility the front, rear and pitch frequency are 1.86 Hz, 2.50 Hz, 2.17 
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5.19: The position of the Pareto-optimal front in the feasible search space in case of three bumps 
excitation at car speed 50 km/h 




Time of travel (s) Time of travel (s) 

Figure 5.20: Results for three bumps excitation at a speed 50 km/h 


Hz respectively. In the maximum transmissibility case the respective frequency values are 1.04 Hz, 1.90 
Hz. and 1.50 Hz. 
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Results 

kfs 

a 

a 

a 

Transmissibility 

Max acceleration 
(m/s2) 

Minimum transmissibility 

4.93 

0.12 

2.35 

0.22 

0.76 

7.8 

Maximum transmissibility 



0.74 

0.27 

1.40 

4.0 


lable 5.12: Optimal values of stiffness and damping parameters for three bumps excitation at speed 50 
km/h 


Results for rounded displacement step excitation at 5 km/h 50 km/h 

In the; rounded displacement step excitation we consider both the car speeds 5 km/h and 50 km/h for 
l)oimce motion. The rounded displacement function as mentioned in chapter 4 has 7 = 50. i.e. the 
excitation is more severe. Henceforth, in all the cases we have considered a cubic damping in addition 
to the linear damping term. The results for the bounce motion are now presented. For the car speed of 
5 km/h the Pareto-optimal front is shown in figure 5.21. The comparision of extreme values is shown 
in the figure 5.22. 



Figure 5.21; The Pareto front showing the optimal values in objective space for rounded displacement 
step at 5 km/h 

Prom table 5.13 it can be noticed that the stiffness values are comparatively lower than those present 
in the previous cases. For minimum transnaissibility value the front, rear and pitch frequencies are 1.91 
Hz, 2.21 Hz, 2.05 Hz respectively. For maximum transmissibility case the respective frequency values 
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Fif^un^ 5.22: Results for rounded displacement step excitation at a speed 5 kin/li 
iire I.IB Hz, 1.77 Hz,. and 1.53 Hz respectively. 


Results 

^fs 


Ks 


Transmissibility 

Max acceleration 
(m/s'-') 

Minimum transmissibility 

3.88 

0.19 

2.46 

0.64 

1.31 

8.83 

Maximum transmissibility 

1.11 

0.18 

2.11 

0.54 

1.58 

6.09 


'lable 5.13: Optimal values of stiffness and damping parameters for rounded displacement step excitation 
at. speed 5 km/h 

III case of rounded step excitation at car speed of 50 km/h the Pareto-front has a wider range of 
t rausmissibility and acceleration values as shown in the figure 5.23. The extreme responses are shown 
in the figure 5.24. The extreme optimal solutions are given in the table 5.14 The front, rear and pitch 
fVecpiency for miiiiuuim transmissibility case are 1.54 Hz, 2.49 Hz, 2.03 Hz respectively. The respective 
vahu's for the case of maximum transmissibility are 0.87 Hz, 1.63 Hz, 1.27 Hz. It is observexi that the 
damping coefficient values tend to be minimum and the various responses can be obtained by varying 
the values of the stiffness parameters with higher stiffness values corresponding to higher acceleration 
and lower displacement or transmissibiUty. 
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Figure 5.23: The Pareto front showing the optimal values in objective space for rounded displacement 
.step at 50 km/h 
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Figure 5.24: Results for rounded displacement step excitation at a speed 50 km/h 

Results for pitching naotion with rounded step displacement at a cm speed of 50 km/h 

Next, we consider another set of objectives i.e., pitching amplitude and pitching acceleration of the 
sprung mass. The excitation function is the rounded displacement step function and the car speed is 
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Results 

^fs 


Ks 

■ 

Transmissibility 

Max acceleration 
(m/s2) 

Minimum transmissibility 

4.90 

0.43 

1.61 

0.10 

1.39 

8.85 

Maximum transmissibility 

2.11 

0.10 

0.51 
1 

0.11 

1.98 

4.22 


Table 5.14: Optimal values of stiffness and damping paxameters for rounded displacement step excitation 
at speed 50 km/h 


50 km/h. The obtained Pareto-front for the pitching motion is shown in the figure 5.25 and the two 
I'xtreme responses are shown in the .figure 5.26. 



Figure 5.25: The Pareto front showing the optimal values in case of pitching motion for rounded 
displacement step at 50 km/h 

Prom table 5.15 it can be made out that for pitching amplitude to be smaller the front suspension 
stiffness should be smaller than the rear suspension stiffness. This observation matches with the theory. 
Also, from the figures showing pitching motion, it is noted that for the faster dying down of pitching 
motion the rear damping co-eflBcient should have a high value. The front, rear and pitch frequencies for 
the minimum transmissibility case are 1.03 Hz, 1.17 Hz, 1.11 Hz respectively. The respective frequency 
values for the maximum transmissibility case are 0.87 Hz, 1.04 Hz, 0.95 Hz. 
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I' 5.26: Results for rounded displacement step excitation in case of pitching motion at a car sp('ed 
50 kiu/li 


Results 

kfs 


jua 

& 

Pitching amplitude 
(rad) 

Max acceleration 
(rad/s'^) 

Minimum pitching amplitude 

0.85 

0.11 

0.93 

1.2 

0.028 

7.56 

Miiximum pitching amplitude 

0.86 

0.10 

0.51 

0.20 

0.049 

2.69 


'I’able 5.15: Optimal values of stiffness and damping parameters for rounded displacement step excitation 
ill case of pitching motion at speed 50 km/h 


Results for rolling motion with a rounded step hole at a car speed of 50 km/h 

fn 3-D model we can consider the rolling motion of the sprung mass as well. So, optimization of the 
rolling motion is done for a special case. In this case we assume that, the right wheels of the car fall in a 
ditch which can be modelled as the negative of rounded displacement step function. The left wheels do 
not. experience any kind of road excitation, and hence a rolling motion of the sprung mass takes plac:e. 
'file results of the optimization of rolling motion axe presented now. The Pareto-optimal front for the 
rolling motion is shown in the figure 5.27. The extreme responses are shown in the figure 5.28. 

From table 5.16 we observe that the stiffness values are hitting their lower bounds. The damping 
coefficient values are the deciding factor in minimizing the rolling motion. Higher damping ctoefficients 
lead to lesser rolling. All the three frequencies are very close to the lower bound, i.e., 0.8 Hz. 


5.3 Conclusion 

In the present study a further step is taken with respect to a systematical improvement of vehicle’s dy- 
namic behaviour. Before presenting the conclusion of the present work, it would be worthwhile to look 
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Maximum rolling amplitude (rad) 


Figure 5.27; The Paxeto front showing the optimal values in case of rolling motion for rounded displace- 
nu'ut step at 50 km/h 




Figure 5.28: Results for rounded displacement step excitation in case of rolling motion at a car speed 
50 km/h 

at the present status of the optimization methods in the industry. This will help us better understand 
the aim of the present study. Three very important points refer to the general opinion of designers 
about the optimization procedures. 


• Optimization algorithms are mathematical tools that cannot take into consideration the physical 
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Results 

kfs 


Ua 

^rs 

c? 

Rolling amplitude 
(rad) 

Max accedfu-ation 
(rad/.s'-^) 

Minimum rolling amplitude 

0.51 

4.55 

0.53 

1.93 

0.034 

2.67 

Maximum rolling amplitude 

0.51 

0.12 

0.52 

0.11 

0.059 

0.72 


I ;i.l)l{' 5. IG; Optimal values of stiffness and damping parameters for rounded displacement st(^) excitation 
in ease' of rolling motion at speed 50 km/h 


or (iiigiiieering aspects of the problem and so, they give us mathematical solutions, that are math- 
(jmatically correct, but most of the time are of no practical interest to the engiiuHir. 


• Physical or engineering aspects of the problem must be taken into consideration by incorporating 
appropriate constraints in the formulation of the optimization problem. 


• Optimization algorithms may be fouled very easily by considering a local minimum. To reach the; 
global optimum a number of runs have to be made from different initial points. This loads to high 
computation cost and time. 

With these points in mind, we proceeded to optimize the ride behaviour of the vehicle. From the 
prcisent study following conclusions can be made. 


• GAs can be easily coupled with numerical integration technique to solve the problems of dynamic 
engineering design. In the present study, the combination of RGA and NSGA-II with adaptive 
st«?p-size control Runge-Kutta method has worked desirably. 


• One of the major problems in optimizing the ride behaviour of a vehicle has been the lack of 
optimization tools which are less time consuming, since the evaluation process involves iterative 
umnerical integration of differential equations of motion. This problem has been effectively ov<;r- 
come by the use of genetic algorithms. 


• Most of the conventional techniques are unable to handle the problem in which the search space 
is highly constrained, since the penalty function methods used distort the objective function. A 
good way of constraint handling in the RGA and NSGA-II algorithms has made the use of a 
number of constraints in the present study possible. 
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• This study to(;usses mainly on the application of multi-objective genetic algorithm, i.e. NSGA-II. 
to the dynamic system design. Since, in most of the practical design problems, dynamic sys- 
t(!uis have to be optimized with respect to several conflicting specifications and several conflicting 
si)eciflc,atiou.s and several diiferent design are acceptable as optimal with respect to the same set 
oi spe(;ification, the present work brings forth an effective optimization strategy for the d(!sigii 
engineer. This strategy can be easily applied to many multi-body dynamics problems. 

Also, a few features of the car suspension design which came forwai’d or were reaffirmed are as 
follows : 


• Two confuting objectives of ride are the minimization of vehicle’s displacement and a<H:el(>ratiun 
respectively. The designer has to strike a balance between these two objectives. A set of optimal 
values provide the designer with a choice and he can select a particular solution which meets his 
needs. 


• It is observed that the softer springs lead to lower natural frequencies and lower accel(;rations but 
the displacement values higher as compared to those of the stifier springs. It is desirable to have 
as soft a suspension as possible, but the stiffness values have to be restricted to a lower bound 
because of rattle space constraints and to improve the handling characteristics of the vehicle. 


• It is observed that the higher stiffness values axe restricted by the upper limit of maximum allowed 
frequency. In case of most cars the natural frequency is kept below 2 Hz. Performance cars on 
which ride is sacrificed for the handling benefits of a car suspension, have natural frequency upto 
2.5 Hz. 


• The amount of damping required is a compromise between that needed due to prevent undue 
persistence of a vibration at natural frequency excited by a single disturbance, which also pre- 
vents the build up of excessive amplitude of a forced vibration due to a series of impulses, and 
the uncomfortable fact that the higher the damping force the greater the disturbance fed into 
the sprung mass by any given road irregularity. We usually obtained low values of damping co- 
efficients. However, too low damping cannot be used as it is insuflicient to control wheel hop 
oscillations which compromise road holding behaviour. 
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• Also, another important feature i.e., non-linear damping has been used in the pres(’nt. study. It 
was found that non-linearity of dampers comes into play only in the case of seven? shock loading. 
So. we have used a rounded displacement step to study the effect of cubic damping and then 
optimization has been performed on suspension with non-linear dampers. It is observed that the 
negative non-linear damping coefficient worsens the displac.ement response!. 

As a finaruot(! one might hypothesize that the ride engineer’s goal should generally b(! to eliminat(; 
all vibrations in a vehicle. Even though this will never be possible in a motor vehicle, it do(!S give 
dirt!ction to the development effort. Yet there are two contrary phenomena that must be dealt with, 
hirst, the elimination of one vibration will always expose another lesser annoyance. Second, in the limit, 
elimination of all vibrations is not desirable as vibrations are the source of road feel considered to be 
esscMitial feedback to the driver of a motor vehicle. 

5.4 Scope for the future work 

'riie present work is one step towards the effective use of multi-objective genetic algorithms to the dy- 
namics design problems. It provides a new strategy in which optimization algorithms efftictively couple 
with numerical integration technique to optimize the ride dynamics of vehicle. However, as suggestetl 
t.li(\ vehicle dynamics is a very complex issue and it can be viewed from different perspe(!tives. Also. 
I )t her fields of engineering design should get encouraged to use multi-objective genetic algorithms as an 
optimization tool. A few suggestions are put forth for the future work related to this topic. 


• Other conflicting parameters like ride handling and ride comfort can be considered. Also, we 
can consider the conflict between pitch and the bounce motion of the vehicle for multi-objective 
optimization. 


• GAs can be used in conjunction with finite element models. This way a problem which combines 
sizing (stiffness and damping properties of elements) and geometry (nodal co-ordinates) optimiza- 
tion can be tackled. 


Another field of interest is the inclusion of kinematic control in whicle handling and the iufluenc.e 
of the tyre characteristics on the vehicle’s behaviour. 



5.5 Closure 


111 this chaptdr the results of optimization of a car suspension system are produced. Two dimensional 
nil idol and three dimensional models have been dealt separately. First the results of single objective 
ii])tiiiiizatioii using real-coded GA are presented and then the optimization results of multi-objective 
optimization results using NSGA-II are presented. Then the conclusion of the whole study is given in 
lu iel’. The <;hai)ter ends with some suggestions for future work related to this topic. 
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